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CUFFHANGER: 
KEEPING THE SCOW 


The 'CUFFHANGER' listings published in this 
magazine and subsequent parts bear absolutely no 
resemblance to, and are in no way associated with, 
the computer game called 'CUFF HANGER' re- 
leased for the Commodore 64 and published by 
New Generation Software Limited. 




the number of lives that Willie’s got left. The 
number 48 is loaded into B and the contents 
of the two registers are added together to give 
the ASCII code for the number that needs to 
be printed. 

The asc routine is then called. This is the 
routine which returns the pointer to the image 
data for the numeral needed in ROM. A is 
then loaded with 41 to set the colour of the 
letter and the print routine is called. This 
actually prints the number of lives Willie has 
left. The tune routine is then called which 
plays the tune. Note that the music routine on 
page 966 must now be re-assembled with an 
origin of 60000 — otherwise it will be 
overwritten by the routine that follows it. 


No game is complete without a 
scoring routine, and even Willie 
won’t be satisfied with just the bits 
of his picnic. Enter these routines 
and play the numbers game 

The scoring routine must also scroll on the 
appropriate screen and print up the score, 
level and lives on the screen. 


The following program calls up the appropri- 
ate screen, prints up the score and number of 
lives Willie has left and plays the tune: 


org 58676 
call Isi 
call scp 
Id hi, 119 
Id a, (57343) 
Id b,48 
add a,b 
call asc 
Id a,41 


call tune 


ON THE STARTING BLOCKS 


The call Isi calls the routine that scrolls on the 
screen. This takes care of which screen is 
required, through its elb routine. 

Then the scp routine is called. This is the 
one that actually prints up the score and has 
not been given to you yet, so don’t run the 
program until you have keyed in the scp 
program given below. HL is loaded with 1 19. 
This is print position of the number of lives. 
The print routine is going to be called and the 
HL register is used to carry the print position 
into it. 

A is loaded with the contents of 57,343. 
This is the memory location set aside to hold 


uiy juowo 


org 58174 


org 58217 


print 


org 60000 


tune 
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HL is incremented to move along to the 
next print position — you’ll notice that the 
score is being printed from the high digit to 
the low. IX is incremented to advance the 
score pointer to the next location containing a 
score digit. The score is stored as decimal 
digits — one in each memory byte— from the 
high digit to the low one, up memory. 

The B counter is popped back off the stack 
and decremented by the djnz instruction. And 
if the counter has not been decremented to 
zero the processor jumps back to print up the 
next digit. 

When it has been round the loop five times, 
printing up the five digits of the score, the 
processor drops out of the loop and returns. 

On the Commodore it is necessary to have a 
routine at this stage which works out where 
positions are on the screen. In previous 
routines you have specified X and Y coordi- 
nates. From those coordinates, the appropri- 
ate screen memory location must be 
calculated. 

The formula to do this is simple. There are 
40 character squares across the screen, so you 
must multiply the Y coordinates by 40 and 
add the X coordinate to it, then add on 1024, 
which is the base address of the screen 
memory, to give the appropriate location. 

Memory location $0352 is used to store the 
X coordinate and $0353 is used to store the Y 


coordinate. And locations $0384 and $0385 
are going to be used to store the resulting 
screen address. These locations are in the 
cassette buffer, but the cassette player is not 
going to be used while the game is being 
played and they are convenient RAM loc- 
ations to use to pass parameters in and out of 
this multiplication routine. 


0RG 

20480 

CLC 


LDA 

$0353 

LDA 

$0384 

STA 

$0384 

ADC 

$0386 

LDA 

#0 

STA 

$0384 

STA 

$0385 

LDA 

$0385 

CLC 


ADC 

$0387 

LDX 

#5 

CLC 


LOOP CLC 


ADC 

#4 

R0L 

$0384 

STA 

$0385 

R0L 

$0385 

CLC 


DEX 


LDA 

$0352 

BNE 

LOOP 

ADC 

$0384 

LDA 

$0353 

STA 

$0384 

STA 

$0386 

LDA 

$0385 

LDA 

#0 

ADC 

#0 

STA 

$0387 

STA 

$0385 

LDX 

#3 

LDA 

$0384 

L00PJ CLC 


STA 

$FB 

R0L 

$0386 

LDA 

$0385 

R0L 

$0387 

STA 

$FC 

DEX 


LDY 

#0 

BNE 

L00PJ 

RTS 



The Y coordinate in $0353 is loaded into the 
accumulator and stored directly in the low 
byte of the result, $0384. Then the high byte, 
$0385, is cleared by storing 0 in it. 

Multiplication on an eight-bit processor is 
not an easy matter. It has to be done by a 
combination of shifts, rotates or repeated 
additions. To multiply a number by 40 it is 
easiest to multiply it by 8 — that is, 2f3 — then 
multiply it by 32 — which is 2|5— and add the 
results together (32 + 8 = 40). Multiplication 
by powers of 2 is a relatively simple matter. 
To multiply by 2, you simply have to shift all 
the bits in a register one place to the left. So to 
multiply by 8 you have to shift the bits to the 
left three places — or shift them to the left one 


This is the scp routine called in the routine 
above that actually prints the score on the 
screen: 


org 58939 

call print 

Id hi, 55 

inc hi 

Id ix, 57337 

inc ix 

Id b,6 

pop be 

push be 

djnz seq 

Id a,(ix + 0) 

ret 

Id b,48 

org 58174 

add a,b 

asc 

call asc 

org 58217 

Id a, 41 

print # 


FIGURING THE FIGURES 


HL is loaded with 55, the screen position the 
leading digit of the score will be printed at. 
The IX index register is loaded with 57,338, 
the address of the first byte which contains 
the score. The register can then be used as a 
pointer. 

B is loaded with 5, the number of digits in 
the score. This counter is saved by pushing it 
onto the stack. 

A is then loaded with the contents of the 
memory location pointed to by the IX 
register — in other words, the first byte con- 
taining the score. 48 is added to it to give the 
ASCII code and the asc routine is called to get 
the image data pointer in ROM. 

Then A is loaded with 41 to set the colour 
and print is called to print out the digit. 


MULTIPLICATION 


SCORING 
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place three times — and to multiply by 32 you 
have to shift them five places to the left, or one 
place five times. 

The problem is dealing with the bits that 
are shifted out of the register. 

Fortunately the 6510 chip has a very useful 
instruction that helps you deal with this. It is 
called a rotate. This is subtly different from. a 
straight shift, as it does not lose the bit from 
the end of the register when a shift is made. 
Nor does it simply fill in the free bit at the 
other end with a 0. A rotate fills the empty bit 
with the contents of the carry flag and resets 
the carry flag with the bit that is shoved out of 
the other end of the register. 

So here the index register X is loaded with 
5 — it is going to count the number of rotates, 
or the number of times that the contents of the 
register are going to be multiplied by 2 — and 
the carry flag is cleared. ROL $0384 rotates the 
contents of memory location $0384, or the Y 
coordinate, one place to the left. Note that 
because the carry flag is clear, a zero is put in 
the least significant bit. So the contents of this 
particular memory location are multiplied by 
2 without any extraneous bits and pieces 
being added in. 

Any bit that has been pushed out of the 
register during the rotate is now in the carry 
flag. So if a rotate is performed on the high 
byte of the result in $0385, any overflow will 
be automatically accounted for. The second 
rotate will take the overflow bit from the carry 
flag and put it into the least significant bit of 
the high byte. 



The high byte will not overflow of course. 
It was set to 0, so you would need to do nine 
rotates before any bit entering the least 
significant bit on the first rotate would be 
shifted out of the most significant bit. 

The X register is then decremented and the 
f BNE instruction loops the processor back until 
the contents of the two bytes have been 
shifted left five times — that is, until the Y 



: 

mm 





80 DATA31/I/I 
90 DATA1 7,1 ,83 
100 DATA99,111,114 
110 DATA1 01,58, 48 
120 DATA48,48,48 
130 DATA48,48,31 
140 DATA1 ,2,76 
150 DATA1 05,1 18,1 01 
160 DATA1 15,58,32 
170 DATA32,32,32 
180 DATA32,53,31 
190 DATA2,4,17 
200 DATA5,1 62,1 63 
210 DATA1 64,1 0,8 
220 DATA8,8,1 65 
230 DATA1 66,1 67,31 
240 DATA1 0,5,1 62 
250 DATA1 63,1 64,10 
260 DATA8,8,8 
270 DATA1 65,1 66,1 67 
280 DATA31,5,4 
290 DATA1 7,3,1 49 
300 DATA1 50,31,0 
310 DATA8,1 49,1 50 
320 FORA% = &1 9E6TO 
&1A2D:READ?A%: 
NEXT 

400 FORPASS = 0TO3 
STEP3 

410 P% = &1 A2E 
420 [OPTPASS 
430 .Heading 


coordinate has been multiplied by 32. 

The Y coordinate in $0353 is then loaded 
into a temporary storage location at $0386 
and the next byte, $0387, is set to 0. The 
contents of these two locations are going to be 
shifted left in exactly the same way as before, 
only this time the loop— and the shift— is only 
three times. So $0386 and $0387 end up 
containing the Y coordinate times 8. 

Now the two results have to be added 
together. So the carry flag is cleared, the 
contents of $0384 are loaded into the ac- 
cumulator and the contents of $0386 are 
added to them. The result is stored back in 
$0384, which is where it is needed. 

Any overflow from this addition will now 
be in the carry flag, so when $0385 is loaded 
into the accumulator and added to the cont- 
ents of $0386, any carry is automatically 
taken into account. 

A further 4 is added into the high byte — 
4 x 256 = 1024. This adds in the base address 
of the screen memory. 

Then the X coordinate in $0352 is loaded 
into the accumulator and added to the low 
byte of the result in $0384. The result is 
stored back in $0384. Any overflow from this 
addition is left in the carry flag. 

The high byte of the result is then loaded 
into the accumulator and has 0 added to it. 
Note, this is an add with carry. So if there is a 
carry from the addition of the low byte with 
the X coordinate it will be added in now. If 
not, the contents of the accumulator remain 
unchanged. Either way, they are stored back 
in $0385. Note that the contents of the high 
byte can’t overflow. For that to happen, the 
resulting address would have to be more than 
65,535 which the highest address permissible 
on the 64K Commodore. 

The calculation is now complete, but for 
added convenience the result is copied from 
$0352 and $0353 into the zero-page memory 
locations $FB and $FC. Later, when the 
result needs to be accessed, say by indirect 
indexed address, the instructions can be 
shaved by a byte because zero page addressing 
is used and speed might be of the essence at 
that point in the game. 


440 LDX#0 
450 .Lbl 

460 LDA&1 9E6,X 
470 JSR&1803 
480 INX 
490 CPX#72 
500 BNELbl 
510 RTS 
520 .PtSL 
530 LDA#31 
540 JSR&FFEE 
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550 LDA # 7 
560 JSR&FFEE 
570 LDA # 1 
580 JSR&FFEE 
590 LDA# 17 
600 JSR&FFEE 
610 LDA# 1 
620 JSR&FFEE 
630 LDX#0 
640 .Lb2 
650 LDA&8A.X 
660 CLC 
670 ADC #48 
680 JSR&FFEE 
690 I NX 
700 CPX#6 
710 BNELb2 
720 LDA #31 
730 JSR&FFEE 
740 LDA #12 
750 JSR&FFEE 
760 LDA #2 
770 JSR&FFEE 
780 LDA&89 
790 CLC 
800 ADC #48 
810 JSR&FFEE 
820 RTS 
830 ] NEXT 


HEADINGS 


The DATA which gives the words ‘score’, 
‘lives’, clouds and seagulls is read into a data 
table from &19E6 to &1A2D. And the ma- 
chine code routine that uses it is assembled 
directly after it. 

The assembly language program starts by 
initializing the index register X to 0 to act as a 
loop counter. The byte of data from the data 
table pointed to by the base address, offset by 
X, is loaded into the accumulator and printed 
at the top of the screen by jumping to the 
user-defined graphic subroutine which is 
located at &1803. 

X is then incremented and compared to 72 
to check whether all the data has been output 
to the screen. If it hasn’t, the processor 
branches back to deal with the next byte. If it 
has, the processor drops out of the loop, hits 
the RTS and returns. 


FIGURING 


Loading A with 31 and calling the routine at 
&FFEE, allows you to position the cursor. 
The 7 and 1 are the coordinates. Loading A 
with 17 and calling &FFEE sets the colour. 
Here colour 1 is used, but don’t forget 
that you’ve redefined all the colours earlier. 

X is initialized to 0 again and the number 
in zero-page memory location &8A, offset by 
X, is loaded into the accumulator. This 
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location and the five following it are used to 
store the score, in decimal digits, from the 
high digit to the low digit, up memory. 

The number 48 is added to the figure to 
give the ASCII, which is output to the screen 
by calling &FFEE. 

X is then incremented and the processor 
loops back to print out the next digit — unless 
it is the case that all six of them have been put 
up on the screen. 

The cursor is then repositioned at 12, 2 and 
the number in &89 — which is where the 
number of lives are stored — is converted into 
ASCII and output to the screen. The pro- 
cessor then returns again. 

To test it with the other machine code in 
the machine CALL&1 B32, putting the level of 
the screen in &83. 


PRSC 


PRSCB 


Ql 


This routine scrolls on the appropriate screen 
and prints up the sun, the lives and the 
score — and makes sure that the numbers that 
are going to be displayed on the screen do not 
bump into each other: 


PRSCA 


DECB 

BNE SCPRD 
JSR 30000 
RTS 
NOP 
NOP 

PSHS D,X,Y 
LDX #18240 
LDB #6 
LDY #1814 
LDA ,X 
PSHS X,B 
BITB #1 
BNE ROLL 
LDB #5 
MUL 

ADDD #17724 
TFR D,X 
PSHS Y 
LDB #5 
LDA ,X + 

STA ,Y 
LEAY 32, Y 
DECB 

BNE PRSCA 
PULS Y 


SCPR 

SCPRI 


SCPRZ 

SCPRC 


THE SCREENING 


The first thing this routine does is jump to the 
subroutine at $4AA5 which scrolls on the 
screen and prints up the sun. 

Then X is loaded with the screen position 
which you want to start printing with the 
word ‘SCORE’. U, the user stack pointer, is 
loaded with the start position of the word data 
in memory. This was input as a data table in 
an earlier part. And B is set to 5. The letters 
are only five bytes deep. 

There are five letters in the word ‘SCORE’ 
plus a space, making six. But A is only loaded 
with 3 as the Y register is going to be used to 
print two bytes at a time. 

The last two bytes of data are pulled off the 
stack — which is, in fact, the area of the data 
table pointed to by U. These two bytes are 
stored at the screen position pointed to by the 
X register and the one next to it. The X 
register is then incremented twice to move it 
onto the next screen position. 

A is then decremented and the BNE SCPRI 
branches back until all six bytes of data, which 
make up a single line of the word, have been 


SCPRD 


ORG 19489 

LEAY 1,Y 

printed. When one line has been printed, 

JSR S4AA5 

ROLRET PULS B,X 

LEAX 26, X increments the X register by 26 to 

LDX #1807 

LEAX 1 ,X 

move it from the right-hand end of one line to 

LDU #17544 

DECB 

the left-hand end of the one below it. There 

LDB #5 

BNE PRSCB 

are 32 locations across the screen and six 

LDA #3 

PULS Y,X,D 

letters: 32 — 6 = 26. 

PULU Y 

RTS 

B is then decremented and the processor 

STY ,X+ + 

ROLL LDB #5 

branches back to deal with the next line if all 5 

DECA 

MUL 

lines of the word have not been printed yet. 

BNE SCPRI 

ADDD #17724 

When they have, JSR PRSC jumps off to the 

LEAX 26, X 

TFR D,X 

subroutine that prints the score on the screen. 

DECB 

PSHS Y 


BNE SCPR 

LDB #5 

■ THAT'S LIVESB * 

JSR PRSC 

RLLA LDA ,X + 

The next routine which prints up the word 

LDX #2063 

PSHS A 

‘LIVES’ is almost an exact repeat of the 

LDU #17574 

ANDA #15 

routine that printed ‘SCORE’ above — only 

LDB #5 

LSLA 

the screen pointer in X and the data pointer in 

LDA #3 

LSLA 

U are loaded with different values. The word 

PULU Y 

LSLA 

is to be printed in a different position, 

STY,X+ + 

LSLA 

obviously, and its data is in a different part of 

DECA 

ORA #5 

the data table. How the printing is done, 

BNE SCPRC 

STA 1,Y 

though, is exactly the same. 

LEAX 26, X 

PULS A 

This routine does not jump to the subrout- 

DECB 

LSRA 

ine to print up the actual number of lives 

BNE SCPRZ 

LSRA 

though. That follows on directly after the end 

LDA 18239 

LSRA 

of the word-print routine. 

LDB #5 

LSRA 

The number of lives Willie has left is 

MUL 

ORA #$50 

stored in 18,239. The contents of this location 

ADDD #17724 

STA ,Y 

are loaded into A. B is loaded with 5 and the 

TFR D,U 

LEAY 32, Y 

contents of the two registers are multiplied 

LDX #2070 

DECB 

together. 

LDB #5 

BNE RLLA 

The display data for the figure has to be 

PULU A 

PULSY 

located in the data table before it can be 

STA ,X 

LEAY 2,Y 

printed on the screen. Each figure requires 

LEAX 32, X 

BRA ROLRET 

five bytes of data to describe it, so to locate the 
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beginning of the right area of data you have to 
count along the data table in multiples of five. 

The result of a MUL — which multiplies the 
contents of A and B — is put in D. 1 7,724 — the 
base address of the table— is then added to 
this, so the result points to the start of the 
appropriate figure’s display data in memory. 

This pointer is then transferred into the 
user stack pointer, U, so that this area of data 
effectively becomes the user stack. 

Again X is loaded with the screen position 
where the number of lives is to be printed and 
B is loaded with 5, to count the five bytes of 
data needed to make up the figure. The first 
byte to figure display data is pulled off the 
user stack into A. Then it is stored at the 
screen position given by X. 

X is incremented by 32 this time, as only 
one figure is to be printed, and the screen 
pointer has to be moved down one screen 
position to fill in the next line of pixels of the 
figure. 

B is decremented and the processor bran- 
ches back to deal with the next byte if all five 
bytes haven’t been printed on the screen yet. 

Then, to finish off, the processor jumps to 
the subroutine at 30,000 which plays the 
tune. And when it returns here, the RTS sends 
it back to BASIC. 

This RTS will be overwritten later, when 
the whole program is put together. It is 
followed by two NOPs — which are No 
OPeration instructions. These do nothing and 
are only here to make enough room for a J SR 
instruction, when the RTS is overwritten. 
Then the processor will have to be directed to 
the action loop driving the whole program. 


THE NUMBERS GAME 


The data for each pixel line of the figures 
takes up a whole byte. So if you print them 
directly next to each other, they tend to merge 
into each other and become illegible. 

To overcome that difficulty, the score print 
routine rolls figures along half a byte to create 
a gap between them and make the score 
readable. 

The PRSC routine begins with the contents 
of D,X,Y being pushed onto the hardware 
stack. This is to preserve their values. In fact, 
it is only the value in Y that is going to be 
needed later. But when you are writing 
machine code it is best to push the contents of 
all the registers you are going to use in a 
subroutine. You might decide that you need 
to carry an important parameter in a register 
later. So the rule is — if in doubt, push it. 

Memory location 18,240 holds the first 
byte of the score. The decimal value of each of 
the six digits of the score — from the highest to 
the lowest — are held in six memory locations 


from 18,240 up memory to 18,246. 

There are six figures to be printed so B is 
loaded with 6. Y is loaded with the print 
position. 

A is loaded with the contents of the 
location pointed to X, which points to the 
score memory. Then the values in X and B are 
preserved by pushing them onto the stack. 

Bit zero of B is then checked. If B is even, 
that is, bit zero is not set, BN E ROLL branches 
off to the ROLL routine, which shoves every 
other figure half a space to the right. But if B 
is odd and bit zero is set, the processor 
continues with the next instruction. 

To print the figure up on the screen, you 
have to proceed in the same way as printing 
the lives. You multiply the figure required by 
5 and add the base address. 

This time, though, the resulting pointer is 
transferred into X and the print position 
pointer in Y is stored by pushing it onto the 
stack. 

The relevant byte of display data for the 
figure in question is then loaded into A and X 
is incremented. And it is then stored at the 
screen position pointed to by Y. Y is in- 
cremented by 32 to move down one line of 
pixels. B is decremented and the processor 
branches back to deal with the next byte of 
data until all five have been output to the 
screen. 

When it has finished the initial screen 
position pointer is pulled off the stack again. 
It is incremented to move it onto the next 
print position to the right. 

The values of counter in B and the score 
memory pointer in X are then pulled off the 
stack again. X is incremented ready to deal 
with the next figure to the right. And B is 
decremented, to count down the figures that 
have been dealt with. The processor then 
branches back to deal with the next figure, if 
they haven’t all been dealt with yet. 

If all the figures have been printed on the 
screen, the Y,X and D registers are restored 
by pulling the values they had at the beginn- 
ing of the subroutine off the stack again. Then 
the processor returns to the place it was called 
in the main routine. 


ROLLING ALONG 


If the figure in B was even and the ROLL 
routine was called, the processor starts off by 
locating the start of the appropriate figure’s 
data in the data table. It multiplies the score 
figure by five and adds on the base address, 
transfers the result into X, pushes the print 
position and loads the B register with the byte 
counter. 

LDA,X 4- loads the accumulator with a byte 
of the figure data and increments the pointer. 



This byte is pushed onto the hardware stack 
to preserve it. 

The contents of A are then ANDed with 15, 
which preserves the low nybble and clears the 
high one. Four LSLAs — Logic Shift Lefts on 
the Accumulator — move the low nybble a bit 
at a time into the high nybble. 

ORing the result with 5 puts yellow— the 
background colour — into the low nybble. 
Then it is stored at the print position given by 
the screen pointer Y plus 1 . This prints the 
right-hand half of the bit pattern on the 
screen at the position two to the right of the 
preceding figure — but shoved half a space to 
the left. 

The whole bit pattern is then pulled off the 
hardware stack again and the high nybble is 
shifted right, into the low nybble. ORing with 
50 hex sets the high nybble to bright yellow. 
And the resulting byte is stored at the print 
position pointed to by Y. This prints the left- 
hand side of the bit pattern in the right-hand 
half of the print position following the pre- 
ceding figure — and marries it up to the right- 
hand half of the bit pattern in the left-hand 
half of the next print position. The other two 
halves of the print positions are filled in with 
the background colour yellow, effectively 
leaving a half digit gap between each figure. 

LEAY 32, Y increments Y by 32 to deal with 
the next line of pixels down the figure. B is 
decremented and the processor branches back 
to pick up the next byte of the bit pattern, 
unless all five lines of pixels have already been 
printed. 

When all that has been done, the processor 
drops out of the loop and increments Y by a 
further 2. This moves the screen pointer two 
positions to the right — effectively the even- 
numbered figure has occupied two adjacent 
screen positions. 

BRA ROLRET then branches always back to 
the label ROLRET, where the pointers are 
updated before the routine jumps back to deal 
with the next — odd-numbered — digit. 
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FOX AND 

A 






Allow your computer to collect its 
thoughts before the game 
commences. With these routines it 
can play fox or geese, and can look 
ahead to improve its play 


This is the third and final part of the Fox and 
Geese article. The remaining routines are 
those which will enable the computer to take 
the part of either the fox or the geese. In fact 
the machine can be set to play both parts, 
when it competes against itself, or you could 
set it to play neither and play against one of 
your friends using the computer instead of a 
board and pieces. 


210 LET P = B(G(1)) + B(G(2)) + B(G(3)) + 
B(G(4)) 

220 LET X = F: IF P<B(32) THEN LET 
P = P — BX: LET X = 33 — F 
230 LET P = P*B(X): RETURN 
250 LET F = FN A(ABS P) — 30: LET 
B = P/B(F): IF B<0 THEN LET 
B = B + BX: LET F = 33 — F 
260 FOR A = 1 TO 4: LET G(A) = FN 
A(B) + 1 : LET B = B - B(G(A)): NEXT A: 
RETURN 

KH 

210P = B(G(1)) + B(G(2)) + B(G(3)) + 
B(G(4)) 

220 X = FrIFP < B(31 )THENP = P — BX: 

X = 31 — F 

230 P = P‘B(X):RETURN 

250 F = FNA(ABS(P)) — 31 :B = P/B(F): 

IFB < 0THEN B = B + BX:F = 31 — F 
260 FORA = 1 T04:G(A) = FNA(B): 

B = B — B(G(A)):NEXTA:RETURN 


210 P = B(G(1)) + B(G(2)) + 

B(G(3)) + B(G(4)) 

220 X = F:IF P< B(31) THEN 
P = P — BX:X = 31 — F 
230 P = P'B(X):RETURN 
250 F = FNA(ABS(P)) — 31 :B = P/B(F):IF 
B<0 THEN B = B + BX:F = 31 — F 
260 F0RA = 1T04:G(A) = FNA(B): B = B — B 
(G(A)):NEXTA:RETURN 


E2 


210 P = B(G(1)) + B(G(2)) + B(G(3)) + 
B(G(4)) 

220 X = F:IFP < B(31 ) THENP = P - BX: 

X = 31 — F 

230 P = P*B(X):RETURN 

250 F = FNA(ABS(P)) — 31 :B = P/B(F): 

IFB < 0 THENB = B + BX:F = 31 — F 
260 FORA = 1T04:G(A) = FNA(B): 

B = B — B(G(A)):NEXT:RETURN 

These are two of the most important routines 
in the program. The routine from Line 210 to 
Line 230 evaluates the playing position and 
packs it into a single number. Conversely, 
when the program has decided on the value of 
the best move, it has to convert that number 
into a move — or set of positions for the pieces. 
Lines 250 and 260 unpack the single value 
into the set of values needed for positioning 
the pieces. 

These subroutines are called very frequ- 
ently during the program, so they have been 
placed right at the start of the program. 



FOX'S MOVE 


1010 GOSUB 210 

1020 GOSUB 310: GOSUB 250 

1030 IF F > 28 THEN PRINT AT 21,0; 

“THE FOX HAS WONDD 

□ □□□□□□□□□□□ 

GOTO 1410 

1032 GOSUB 410: IF V = H THEN PRINT 
“THE GEESE HAVE WONDD 

DDDDDDDDDDD”: 

GOTO 1410 

1040 IF PF THEN GOTO 1110 
1050 INPUT “MOVE FOX TO ”;B: IF B= — 1 
THEN GOSUB 2710: GOTO 1030 



X = B THEN IF NOT (FN X(X)) THEN LET 
F = B: GOSUB 210: 

GOTO 1200 

1070 NEXT A: PRINT AT 21,0;“THAT 
IS NOT LEGALDDDDD 
DDDDDDDDDD”: GOTO 1050 
1110 LET L = SF: LET M = SF: LET 
V(M) = E*M: IF M>4 THEN DIM 
R(HF + 3): DIM S(HF + 3) 

1112 GOSUB 1120: GOTO 1200 
1120 IF L=1 THEN GOTO 410 
1122 IF L<M — 2 THEN GOSUB 1610: IF 
V< >0 THEN RETURN 
1130 LET L = L — 1 : LET V(L) = H*L: LET 
A(L) = X(F): LET F(L) = F 
1140 LET F = M(A(L),F(L)) 

1150 IF FN X(F) = 0 THEN GOSUB 1320: IF 


1152 
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■ 

PLAYING AT HIGHER LEVELS 

FOX’S MOVE ROUTINE 

■ 

■ 

PLAYING AT HIGHER LEVELS 

USING THE ALPHA-BETA 

■ 

MOVING THE GEESE 


ALGORITHM 

■ 

CHOOSING THE BEST MOVE 

IN A ‘ONE MOVER’ 

■ 

■ 

USING THE HASH-CODE TABLES 

FOR A FASTER GAME 


V>V(L) THEN LET V(L) = V: LET P(L) = F: 
IF V>V(L + 1) THEN LET F = F(L): LET 
L = L + 1: RETURN 

1160 LET A(L) = A(L) — 1: IF A(L) >0 THEN 
GOTO 1140 

1170 LET V = V(L): LET F = F(L): LET 
L = L + 1: IF L = M THEN LET 
F = P(M — 1): GOSUB 210: RETURN 
1172 IF L<M — 2 THEN GOSUB 1510: 

RETURN 
1180 RETURN 

SE 

1010 GOSUB 210 

1020 GOSUB 31 0:GOSUB 250 

1030 IFF > 27THEN PRINTTAB(8); 

“THE FOX HAS WON”:GOTO1410 



1032 GOSUB41 0:IFV = HTHEN PRINT 
TAB(8);“THE GEESE HAVE 
WON”:GOTO1410 
1040 IF PF THEN PRINTTAB(8); 

“THINKING. . ”:GOTO 1110 
1050 INPUT UUUIIUUUII 
MOVE FOX TO”;B 
1055 GOSUB 340:IF B= — 1 THEN 
PRINTTAB(8);“THE GEESE HAVE 
WON”:GOTO1410 

1060 FORA = 0TOX(F):X = M(A,F):IFX = B 
THENIFNOTFNX(X)THENF = B:GOSUB210: 
GOTO1200 

1070 NEXTA:PRINTTAB(8);“THAT IS NOT 
LEGAL”:GOTO1050 

1110 L = SF:M =SF:V(M) = E'M:IFM >4 
THENFORA = 0TOHF:R(A) = 0: 

NEXTA 

1112 GOSUB1 1 20:GOTO1 200 
1120 IFL = 1THEN410 



1122 IFL < M — 2THENGOSUB1 61 0: 

I FV < >0THEN RETURN 
1130 L= L— 1 :V(L) = H‘L:A(L) = X(F): 

F(L) = F 

1140 F = M(A(L),F(L)) 

1150 IFFNX(F)o0THEN 1160 
1 1 55 GOSUB1 320:IFV > V(L)THENV(L) = V: 
P(L) = F:IFV> V(L + 1)THENF = F(L): 

L = L + 1:RETURN 

1160 A(L) = A(L) — 1 :IFA(L) > = 0THEN1 1 40 
1170 V = V(L):F = F(L):L = L + 1:IFL=M 
THENF = P(M - 1 ):GOSUB210:RETURN 
1172 IFL< M — 2THENGOSUB1 51 0:RETURN 
1180 RETURN 


THEN FOR A = 0TOHF:R(A) = 0:NEXTA 
1112 GOSUB 1120:GOTO 1200 
1120 IFL = 1 THEN 410 
1122 IF L<(M — 2) THEN GOSUB 1610: 

IFV< >0 THEN RETURN 
1130 L = L— 1:V(L) = H’L:A(L) = X(F): 

F(L) = F 

1140 F = M(A(L),F(L)) 

1150 IF FNX(F) = 0 THEN GOSUB 1320:IF 
V > V(L) THEN V(L) = V:P(L) = F: 

IFV>V(L + 1) THEN F = F(L):L = L + 1: 
RETURN 

1160 A(L) = A(L) — 1 :IFA(L) > = 0 THEN 1140 
1170 V=V(L):F = F(L):L = L + 1: 

IFL = M THEN F = P(M — 1): 


a 

1010 GOSUB 210 
1020 GOSUB 31 0:GOSUB 250 
1030 IF F>27 THEN PRINT TAB(8);“THE FOX 
HAS WON”:GOTO1410 
1032 GOSUB410:IF V = H THEN PRINT 
TAB(8);“THE GEESE HAVE 
WON”:GOTO1410 
1040 IF PF THEN PRINT TAB(8); 

“THINKING... ”:GOTO 1110 
1050 INPUT TAB(8)“MOVE FOX TOD” B 
1055 GOSUB 340 

1060 A = 0:REPEATX = M(A,F):IFX< >B OR 
FNX(X) THEN A=A + 1:UNTIL 
A > X(F):ELSE UNTIL TRUE: 

F = B:GOSUB210:GOTO1200 
1070 PRINT TAB(8);“THAT’S NOT LEGAL”: 
GOTO 1050 

1110 L = SF:M = SF:V(M) = E’M:IFM>4 


GOSUB210:RETURN 
1172 IF L< (M — 2) THEN GOSUB1510: 

RETURN 
1180 RETURN 




1010 SCREEN1,0:GOSUB210:GOTO1030 
1020 XX = FNXX(G):YY = FNYY(G):PUT 
(XX, YY) - (XX + 19,YY + 19),SQ,PSET:XX 
= FNXX(G(C)):YY = FNYY(G(C)) + 5: 
PUT(XX,YY) - (XX + 1 9,YY + 9),GS,PSET 
1030 CLS:IF F > 27 THEN PLAYV$: 
PRINT@7,“THE FOX HAS 
WON”:GOTO1410 

1032 GOSUB41 0:1F V = H THENPLAYV$: 
PRINT@6,“THE GEESE HAVE 
WON”:GOTO1410 

1040 LINE(180,0) — (255, 191), PRESET, 

BF:IF PF THENDRAW“BM180,50C4” + 
TH$:G0T01 110 
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1050 DRAW“BM180,80C4” + MW$:XX = 
FNXX(F):YY = FNYY(F):G0SUB1 81 0: 

B = 4'INT (YY/20): B = FNCN(B) 

1055 IF B= — 1 TH EN PLAYV$: PR I NT@6, 
“THE GEESE HAVE WON”:GOTO1410 
1060 FORA = 0TOX(F):X = M(A,F): 

IFX = B AND NOTFNX(X) THENA = X(F): 
XX = FNXX( F) : YY = FN YY( F) : PUT 
(XX, YY) — (XX + 1 9,YY + 1 9), 
SQ,PSET:F= B:GOSUB210:NEXT: 


after the move. The number of plies that the 
program is looking at, M, and the current ply 
being considered, L, are set up in Line 1110. 
Lines 1120 to 1180 are a subroutine to 
evaluate the best move. 

Four arrays are used in the best move 
subroutine: A contains the number of moves 
still to try; V, the best result so far; P, the move 
that yields that result; and F, the fox’s 
previous position. 


GOTO1200 

1070 NEXT:GOSUB5000:GOTO1 040 
1110 L = SF:M = SF:V(M) = E'M:IFM>4 
THENFORA = 0TOHF:R(A) == 0:NEXT 
1112 LF=F:GOSUB1120:XX= FNXX(LF): 

YY = FNYY(LF): PUT(XX,YY) - (XX + 1 9, 

YY + 1 9),SQ,PSET :G0T01 200 
1120 I FL = 1 THEN410 
1122 IFL < M — 2 G0SUB1 61 0:IF Mod) 
THEN RETURN 

1130 L=L — 1:V(L) = H*L:A(L) =X(F): 

F(L) = F 

1140 F = M(A(L),F(L)) 

1 1 50 IFFNX(F) = 0 G0SUB1 320:1 FV > V(L) 
THENV(L) = V:P(L) = F:IF V>V(L + 1) 

THEN F = F(L):L = L + 1: 

RETURN 

1160 A(L) = A(L) — 1 :IF A(L)> = 0THEN 
1140 

1170 V = V(L):F = F(L):L = L + 1: 

IFL=M THENF = P(M — 1): 
GOSUB210:RETURN 
1172 IF L<M — 2 GOSUB 1510 
1180 RETURN 

Lines 1020 to 1180 handle the fox’s move. 
The routine uses the board display subrout- 
ine to display the current status of the board, 
then goes on to check if the fox has won in 
Line 1030. It then checks if there is at least 
one legal move open to the fox (or else the 
geese win) in Line 1032. 

If the player is controlling the geese, Lines 
1050 to 1070 take the input and check its 
legality. If, on the other hand, the computer is 
controlling the fox, Lines 1110 to 1112 look 


MOVING THE GEESE 


500 LET V = 0: FOR C = 1 TO 4: LET 
G = G(C): FOR A = 1 TO Z(G): LET 
X = M(A,G): IF X< > F THEN IF NOT FN 
X(X) THEN RETURN 
502 NEXT A: NEXT C: LET V = 1 : 

RETURN 

1200 GOSUB 310: GOSUB 250 
1202 IF F > 28 THEN PRINT AT 21,0; 

“THE FOX HAS WONDDD 

□ □□□□□□□□□□□□”: 

GOTO 1410 

1204 GOSUB 500: IF V THEN PRINT AT 
21 ,0;“THE FOX HAS WON □ □ □ □ □ 

□ □□□□□□□□□□”: GOTO 
1410 

1210 IF PG THEN GOTO 1310 
1220 INPUT “WHICH GOOSE TO MOVE ? ”;G: 
IF G= -1 THEN GOSUB 2710: GOTO 
1202 

1230 LET C = FN Z(G): IF C = 0 THEN PRINT 
AT 21 ,0;“NO GOOSE AT ’’;G;“”: GOTO 
1220 

1240 INPUT “WHERE TO ”;l: IF 1= -1 THEN 
GOSUB 2710: GOTO 1202 
1250 IF FN X(l) THEN PRINT AT 21,0; 

“NOT ONTO ANOTHER GOOSE 

□ □□□□□□□”: GOTO 1220 
1260 IF l = F THEN PRINT AT 21,0; 

“NOT ONTO FOXD □□□□□□ 

□ □□□□□□□□□□ 

GOTO 1220 


1270 FOR A = 1 TO Z(G): IF M(A,G) = I THEN 
LET G(C) = I: GOTO 1010 
1280 NEXT A: PRINT AT 21 ,0;“THAT 
IS NOT LEGALD □□□□□□ 

□ □□□□□□□”: GOTO 1220 
1310 LET L=SG: LET M = SG: LET 
V(M) = H'M: IF M >4 THEN DIM 
R(HF + 3): DIM S(HF + 3) 

1312 GOSUB 1320: GOTO 1020 
1320 IF L = 1 THEN GOTO 510 
1322 IF L<M — 2 THEN GOSUB 1610: IF 
V< >0 THEN RETURN 
1324 LET L=L — 1: LET V(L) = E*L: LET C = 1 
1330 LET C(L) = C: LET F(L) = G(C): LET 
A(L) = 1: IF A(L)>Z(G(C)) THEN GOTO 
1362 

1340 LET B = M(A(L),F(L)): LET X = FN X(B): 
LET G(C) = B: IF X OR B = F THEN GOTO 
1360 

1350 GOSUB 1120: LET C = C(L): IF V< V(L) 
THEN LET V(L)=V: LET 
P(L) = G(C) + C"32 

1355 IF V < V(L) THEN LET G(C) = F(L): LET 
L = L + 1: RETURN 



A(L) < =Z(F(L)) THEN GOTO 1340 
1362 LET G(C) = F(L): LET C = C + 1: IF 
C<5 THEN GOTO 1330 
1370 LET V = V(L): LET L=L + 1: IF L=M 
THEN LET C = INT (P(L — 1 )/32): LET 
G(C) = P(L — 1 ) — C*32: GOSUB 210: 
RETURN 

1372 IF L<M — 2 THEN GOSUB 1510: 

RETURN 
1380 RETURN 

KH 

500 V = - 1 :FORC = 1T04:G = G(C):IF 
Z(G) <0THENNEXT:RETURN 
502 FORA = 0TOZ(G):X = M(A,G):V = VAND 
(FNX(X)ORX = F):NEXT:NEXT:RETURN 
1200 GOSUB31 0:GOSUB250 
1202 IFF > 27TH EN PRINTTAB(8); 

“THE FOX HAS WON”: 

GOTO1410 

1204 GOSUB500:IFVTHEN PRINT 
TAB(8);“THE FOX HAS WON”: 
GOTO1410 


1154 









1210 IF PG THEN PRINTTAB(8); 

“THINKING... ”:GOTO 1310 
1220 INPUT“UUUUUfeJUU 

WHICH GOOSE TO MOVE”;G 
1225 GOSUB 340 

1 230 C = FNZ(G):IFC = 0THENPRINT 
TAB(8);“NO GOOSE AT”;G: 
GOTO1220 

1240 INPUT“UUUIIIIUIIU 

WHERE TO”;l 

1245 GOSUB 340:IF I =-1 THEN 
PRINTTAB(8);“THE FOX HAS 
WON”:GOTO1410 
1250 IFFNX(I)THEN PRINTTAB(8); 
“NOT ONTO ANOTHER GOOSE”: 
GOTO1220 

1260 IFI = FTHEN PRINTTAB(8); 

“NOT ONTO FOX”: 

GOTO1220 


1270 FORA = — 1TOZ(G):IFA> = 0THENIF 
M(A,G) = ITHENG(C) = 1: 

GOTO1010 



1280 NEXTA:PRINT TAB(8);“THAT IS NOT 
LEGAL”:GOTO1220 

1310 L = SG:M = SG:V(M) = H‘M:IFM >4 
THEN FORA = 0TOHF: 

R(A) = 0:NEXTA 
1312 GOSUB1320:GOTO1 020 
1320 IFL = 1THEN51 0 

1322 IFL < M — 2THENGOSUB1 61 0:IFV < >0 
THENRETURN 

1324 1= L— 1:V(L) = E’L:C = 1 
1330 C(L) = C:F(L) = G(C):A(L) = 0:IF 
A(L) >Z(G(C))THEN1362 
1340 B = M(A(L),F(L)):X = FNX(B):G(C) = B: 

IFXORB = FTHEN1 360 
1350 GOSUB1 1 20:C = C(L) 

1355 IFV<V(L)THENV(L) =V:P(L) = G(C) + 
C‘32:IFV < V(L + l)THENG(C) = F(L): • 

L = L + 1:RETURN 

1 360 A(L) = A(L) + 1 :IFA(L) < = Z(F(L))THEN 
1340 

1 362 G(C) = F(L):C = C + 1 :IFC < 5THEN1 330 
1370 V = V(L):L=L + 1:IFL=MTHENC = INT 


(P(L — 1 )/32):G(C) = P(L — 1 ) AND31 : 
GOSUB210:RETURN 
1372 IFL < M — 2THENGOSUB1 51 0: 

RETURN 
1380 RETURN 

LI 

500 V= — 1:FOR C = 1T04:G = G(C):IF 


1 350 GOSUB1 1 20:C = C(L):IFV < V(L)THEN 
V(L) = V:P(L) = G(C) + C*32:IF 
V<V(L + 1) THEN G(C) = F(L):L = L + 1: 
RETURN 

1360 A(L)=A(L) + 1:IFA(L)< =Z(F(L)) 
THEN1340 

1 362 G(C) = F(L):C = C + 1 :IF C < 5 THEN 
1330 



500 V= — 1 : FORC = 1 T04:G = G(C):IF 
Z(G) <0THENNEXT:RETURN 
502 FORA = 0TOZ(G):X = M(A,G):V = V AND 
(FNX(X)ORX = F):NEXTA,C:RETURN 
1200 GOSUB250:XX = FNXX(F): 

YY = FNYY(F) + 5:PUT(XX,YY) - (XX + 1 9, 
YY + 8),FX,PSET 

1202 IF F>27THENPLAYV$:PRINT 
@6, “THE FOX HAS WON”: 

GOTO1410 

1204 GOSUB 500:IFV THEN PLAYV$: 
PRINT“THE FOX HAS WON”: 

GOTO1410 

1210 LINE(1 80,0) - (255,1 91 ), PRESET, BF:I F 
PG THENDRAW“BM180,50C4” + TH$: 
GOTO1310 

1220 XX = FNXX(G(1 )):YY = FNYY(G(1 )): 

DRAW“BM1 80,80C2” + WG$:GOSUB1 81 0 
1225 G = 4*INT(YY/20):G = FNCN(G) 

1 230 C = FNZ(G):IF C = 0 GOSUB5000: 
GOTO1210 

1 240 DRAW“BM1 80,1 1 0C3” + MW$:GOSUB 
1810:1 = 4’INT(YY/20):I = FNCN(I) 

1245 IF 1 = 1 THENPLAYV$:PRINT@7, 

“THE FOX HAS WON”:GOTO1410 
1250 IFFNX(I) GOSUB5000:GOTO1210 
1260 IFI = F GOSUB5000:GOTO1 21 0 
1270 FORA= — 1TOZ(G):IFA> =0 THENIF 
M(A,G) = I THENXX = FNXX(G(C)): 

YY = FNYY(G(C)):PUT(XX,YY) — (XX " 


Z(G) <0 THEN NEXTC:RETURN 
502 FOR A = 0 TO Z(G):X = M(A,G): 

V = V AND (FNX(X) OR X = F):NEXTA: 
NEXTC:RETURN 
1200 GOSUB310:GOSUB250 
1202 IF F>27 THEN PRINT TAB(8); 

“THE FOX HAS WON”:GOTO1410 
1204 GOSUB 500: IF V THEN PRINT 

TAB(8);“THE FOX HAS WON”:GOTO1410 
1210 IF PG THEN PRINT TAB(8); 

“THINKING... ”:G0T01 310 
1220 INPUTTAB(8)“WHICH GOOSE TO 
MOVED” G:IF G = — 1 THEN 
GOSUB2710: GOTO1202 
1230 C = FNZ(G):IFC = 0 THEN PRINT 
TAB(8);“NO GOOSE ATD”;G:GOTO 1220 
1240 INPUT TAB(8)“WHERE TO”, I 
1250 IFFNX(I) THEN PRINTTAB(8); 

“NOT ONTO ANOTHER GOOSE”: 

GOTO 1220 

1260 IFI = FTHEN PRINT TAB(8); 

“NOT ONTO FOX”:GOTO1220 
1270 IF Z(G)>=0 THEN A = 0: 

REPEAT :IF M(A,G) < > I THEN 
A = A + 1:UNTILA>Z(G):ELSE UNTIL 
TRUE:G(C) = I:GOTO1010 
1280 PRINTTAB(8);“THAT’S NOT 
LEGAL”:GOTO1220 
1310 L = SG:M = SG:V(M) = H'M: 

IFM >4 THEN FOR A = 0 TO HF: 

R(A) = 0:NEXTA 
1312 GOSUB1 320:GOTO1 020 
1320 IFL = 1 THEN510 
1322 IF L< M — 2 THEN GOSUB1610: 

IFV< >0 THEN RETURN 
1324 L = L— 1:V,(L) = E*L:C = 1 
1330 C(L) = C:F(L) = G(C):A(L) = 0:IF 
A(L) > Z(G(C))THEN1 362 
1340 B = M(A(L),F(L)):X = FNX(B): 

G(C) = B:IF X OR B = FTHEN 1360 


1370 V = V(L):L=L + 1:IFL = M THEN 
C = INT(P(L — 1 )/32):G(C) = P(L — 1 ) AND 
31 :GOSUB21 0:RETURN 
1372 IFL < M — 2 THEN GOSUB1510: 

RETURN 
1380 RETURN 
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YY + 1 9),SQ,PSET:G(C) = l:XX = FNXX(I): 
YY = FNYY(I):PUT(XX,YY + 5) - (XX + 1 9, 
YY + 14),GS,PSET:A = Z(G):NEXT: 
GOTO1010 

1280 NEXT:GOSUB5000:GOTO1210 
1310 L = SG:M = SG:V(M) = H*M:IFM>4 
THENFORA = 0TOHF:R(A)=0:NEXT 
1312 GOSUB1320:GOTO1 020 
1320 IFL = 1 THEN510 
1322 IFL < M — 2 GOSUB1610:IF V< >0 
THEN RETURN 

1324 L = L— 1:V(L) = E*L:C = 1 
1330 C(L) = C:F(L) = G(C):A(L) = 0:IF 
A(L)>Z(G(C)) THEN1362 
1340 B = M(A(L),F(L)):X = FNX(B):G = G(C): 

G(C) = B:IF X ORB = F GOTO1360 
1350 G0SUB1 1 20:C = C(L):1FV < V(L)THEN 
V(L) = V:P(L) = G(C) + C*32:IF 
V<V(L + 1)THENG = G(C):G(C) = F(L): 
L=L+1:RETURN 

1360 A(L) = A(L) + 1:IFA(L) < =Z(F(L)) 
THEN1340 

1362 G = G(C):G(C) = F(L):C = C + 1:IFC<5 
THEN1330 

1 370 V = V(L):L = L + 1 :IFL = M THENC = INT 
(P(L— 1 )/32):G = G(C):G(C) = P(L- 1 ) 
AND31:GOSUB210:RETURN 
1372 IFL < M — 2 G0SUB1 51 0 
1380 RETURN 

The routine from Line 1200 to Line 1380 
handles the geese. There are subroutines for 
when the player controls the geese — lines 
1220 to 1290— and when the computer con- 
trols the geese — Lines 1320 to 1380. 




410 LET V = H: FOR A = X(F) TO 1 STEP -1: 
LET X=M(A,F): IF FN X(X) THEN NEXT A: 
LET L = 1 : RETURN 

420 LET B = F: LET F = X: GOSUB 210: LET 

V = P: LET F = B: LET L = 1 : RETURN 
510 LET V = E: FOR C = 1 TO 4: LET 

G = G(C): IF — B(G) > V THEN GOTO 530 
520 FOR A = 1 TO Z(G): LET X = M(A,G): IF 
FN X(X) OR (X = F) THEN NEXT A: GOTO 
530 

528 LET V = B(X) — B(G): LET D = C: LET 

B _ x 

530 NEXT C: LET G = G(D): LET G(D) = B: 
GOSUB 210: LETV = P: LET G(D) = G: 
RETURN 

SB 

410 V = H:FORA = X(F)TO0STEP — 1: 

X = M(A,F):IFFNX(X) < 0THENNEXT:L = 1 : 
RETURN 

420 B = F:F = X:GOSUB21 0:V = P:F = B: 

L = 1:RETURN 

510 V = E:F0RC = 1T04:G = G(C):IF 
— B(G) > VTHEN530 

520 FORA = 0TOZ(G) :X = M(A,G):I FFNX(X) 
ORX = FTHENNEXT:GOTO530 
528 V = B(X) — B(G):D = C:B = X 
530 NEXTC:G = G(D):G(D) = B:GOSUB210: 

V = P:G(D) = G:RETURN 


410 V = H:A = X(F):REPEAT:X= M(A,F):IF 
FNX(X) <0 THEN A = A — 1:UNTIL A<0: 
L=1:RETURN 

420 UNTILTRUE:B = F:F = X:GOSUB 210: 

V = P:F=B:L = 1:RETURN 
510 V=E:FOR C = 1T04:G = G(C):IF 
— B(G)>V OR Z(G) <0THEN 530 
520 A = 0:REPEAT:X = M(A,G):IF FNX(X) OR 
(X = F) THEN A = A + 1:UNTIL A > Z(G): 
GOTO530 


41 0 V = H : FORA = X( F)TO0 STEP - 1 : 

X = M(A,F):IFFNX(X) < 0 THENNEXT:L = 1 : 
RETURN 

420 B = F:F=X:GOSUB210:V = P:F = B: 

L = 1:A = 0:NEXT:RETURN 
510 V = E:F0RC = 1T04:G = G(C):IF 
— B(G) > V THEN530 

520 FORA = 0TOZ(G):X= M(A,G):IFFNX(X) 

OR X=FTHENNEXT:GOTO530 
528 V = B(X) — B(G):D = C:B = X 
530 NEXT:G = G(D):G(D) = B:GOSUB210: 
V=P:G(D) = G:IFSG = 1 THENG(D) = B: 

C = D 

540 RETURN 

These subroutines are concerned with the 
‘one mover’ — in other words, at this stage the 
computer is only looking one move ahead 
when searching for the best move. 

Lines 410 and 420 go through all the 
possible moves open to the fox, using the map 
array M, set up in the subroutine starting at 
Line 2110 (see earlier). The subroutine re- 
turns a value of P, configuration after best 
move, and V, evaluation after best move. 

Lines 510 and 530 are a similar routine, 
but this time it looks for the best move for the 
geese. P and V are set in the same way as in the 
previous subroutine. 

In either case, GOSUB 210 picks the best 
move from those evaluated. 


THE HASH CODE TABLE 


1510 GOSUB 210: LET C = P 

1520 LET C = C — INT ((C/HF + C) — C)'HF: 

IF C<0 OR C> =HF THEN GOTO 1520 
1550 FOR A = C + 1 TO C+4: IF R(A)< >0 
AND R(A) < > P THEN NEXT A: RETURN 
1560 LET R(A) = P: LET S(A) =V: RETURN 
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1510 GOSUB210:C = P 

1520 C = C — INT((C/HF + C) — C)'HF:IFC<0 
ORC> = HF THEN1520 
1550 FORA = C TOC + C:IFR(A)o0AND 
R(A) < > PTHENNEXT:RETURN 
1560 R(A) = P:S(A) = V:A = C + C:NEXT: 
RETURN 

1610 GOSUB21 0:C = P 

1620 C = C - INT((C/HF + C) - C)*HF:IFC < 0 
ORC> =HF THEN1620 
1650 FORA = C TOC + C:IFR(A) < > 0 AND 
R(A) < > P THENNEXT:V = 0: RETURN 
1660 V= -S(A)*(R(A) = P):A = C + C:NEXT: 
RETURN 


1610 GOSUB 210: LET C = P 

1620 LET C = C — INT ((C/HF+ C) — C)*HF: 

IF C<0 OR C> =HF THEN GOTO 1620 
1650 FOR A = C + 1 TOC + 4: IF R(A) < >0 
AND R(A) < > P THEN NEXT A: LET 
V = 0: RETURN 

1660 LET V=S(A)'(R(A) = P): RETURN 


sa 

1510 GOSUB21 0:C = P 

1520 C = C — INT((C/HF + C) — C)'HF:IFC<0 
ORC> = HFTHEN1520 
1 550 FORA = CTOC + 4:IFR(A) < > 0AND 
R(A) < > PTHENNEXT:RETURN 
1560 R(A) = P:S(A)=V: RETURN 
1610 GOSUB210:C = P 
1620 


ORC> = HFTHEN1 620 
1650 FORA = CTOC + 4:IFR(A) < >0AND 
R(A)< >PTHENNEXT:V = 0: RETURN 
1660 V= — S(A)*(R(A) = P): RETURN 
2500 DIM(1999), S(1 999) 

II 

1510 GOSUB 21 0:C = P 
1512 IF ABS(C/HF) >1E9 THEN 
C = C — ((C/HF + C) — C)'HF: 
G0T01512 

1520 C = C-INT((C/HF + C)-C)*HF: 

IF C<0 OR C> = HF THEN 1520 
1550 A = C:REPEAT:IF R(A)< >0 AND 
R(A)< >P THEN A = A + 1:UNTIL 
A = C + 4:RETURN:ELSE UNTIL TRUE 
1560 R(A) = P:S(A) = V:RETURN 
1610 GOSUB 21 0:C = P 
1612 IF ABS(C/HF) >1E9 THEN 
C = C — ((C/HF + C) — C)*HF: 
G0T01612 

1620 C = C-INT((C/HF + C)-C)*HF: 

IF C<0 OR C> =HF THEN 1620 
1650 A = C:REPEAT:IF R(A(< >0 AND 
R(A)< >P THEN A = A + 1:UNTIL 
A = C + 4:V = 0:RETURN:ELSE UNTIL 
TRUE 

1660 V= — S(A)*(R(A) = P):RETURN 


1 81 0 SCREEN1 ,0: PUT(XX,YY) - (XX + 1 9, 

YY + 1 9),SQ, NOT: FORZ = 1 TOI 00: N EXT 
1 820 PUT(XX,YY) - (XX + 1 9,YY + 1 9),SQ, 

NOT 

1830 K$ = INKEY$:IF K$ = “T” AND YY>8 
AND XX >8 THEN YY = YY-20:XX= 

XX — 20:GOTO1810 

1840 IF K$ = CHR$(10) AND YY < 1 29 AND 
XX <129 THEN YY = YY + 20: 

XX = XX + 20:GOTO1 81 0 
1850 IF K$ = CHR$(8) AND XX >28 THEN 
XX = XX — 40:GOTO1 81 0 
1860 IF K$ = CHR$(9) AND XX <128 
THENXX = XX + 40:GOTO1810 
1870 IF K$ = CHR$(13) THEN RETURN 
1875 IF K$ = “Q” THEN YY = 0:XX = -12: 
RETURN 

1880 GOTO 1810 

In the higher levels of play, you’ll want to 
apply the alpha-beta algorithm — see page 
1098. In fact, you’ve already entered the 
algorithm as part of Fox’s move and Geese 
move routines. Before the algorithm is ap- 
plied, they check if it is appropriate to use the 
algorithm — is the chosen level of play suffi- 
ciently high to warrant its use? 

The routine from Line 1 1 10 to Line 1150 


looks after the fox, while the routine from 
Line 1310 to 1350 looks after the geese. 

The algorithm is applied in the last IF test 
at the end of Lines 1 1 50 and 1 350, after V( M ) 
has been set to the appropriate level in Lines 
11 10 and 1310. 

The alpha-beta algorithm is most efficient 
if the computer considers what are likely to be 
the best moves first — for the geese, the 
highest-numbered square in each row of four 
squares, and for the fox, the square open to it 
that is nearest the geese end. 

The alpha-beta algorithm is used in con- 
junction with a technique known as hash- 
coding to build up and use a table of the moves 
already considered. Hash-coding allows the 
computer to check quickly if the move has 
already been considered. The larger the hash- 
code table that can be built into the program, 
the faster-running that program should be. 

The table is initialized in Lines 25#, 2750 
and 2800. There are theoretical ‘best values’ 
for the dimensions of the arrays holding the 
tables — Line 25#. The arrays have been 
DIMensioned to as large as possible, given the 
available RAM in each of the machines. 

The table is zeroed in Lines 1110 and 
1310, the contents are checked in Lines 1122 
and 1322, and set in Lines 1172 and 1372. 
The checking subroutine starts at Line 1610, 
and the setting subroutine starts at Line 1510. 
BBC owners with disk filing systems should 
enter the following before RUNning. 



‘TAPE 

FOR A% = 0TO&1 600:?(&E00 + A%) = ? 

(PAGE + A%):NEXT 
PAGE = &E00 
OLD 
RUN 

The Dragon/Tandy game uses a flashing 
cursor to move the pieces. Mov e it to t he 
square or piece you wish and press |RETURN| in 
response to the prompts to the right of the 
screen. 





PREDICTING TH 



What are the chances of your 
selection scoring a big win on the 
football Pools this weekend? You 
might not improve your luck, but 
may understand why you’ve lost 


Instructing a novice at the controls of a 
modern jetliner or a fighter-bomber would be 
wastefully expensive, when factors such as 
fuel, landing fees and back-up facilities are 
taken into account. This is why military 
authorities find it cheaper to expend large 
sums for trainer aircraft and simulators — 
computer controlled mock-ups that never 
leave the ground, but give a trainee a realistic 
impression of flying. 

The same principle holds true in a large 
number of circumstances — including games 
programming, industrial development and 
marketing, scientific research and govern- 
ment. It is better to program a computer to 
predict the probable result of structural fai- 
lure, or of a particular economic policy, for 
example, than to do the actual experiment, 
which might last five years or even more in 
some circumstances. Not surprisingly, simul- 
ations are a useful and favoured facility. But 
their complexity means that they were only 
really made possible by the advent of com- 
puters. Now, even the home micro can make 
light work of simulation. 

One of the chief tools of such computer 
modelling is the set of rules provided by the 
mathematical study of probability. In this 
context, you may find it helpful to look back 
to the article on pages 694 to 700, which 
explains some of the theory behind this. The 
rest of the model comes from an analysis of 
statistical information, gathered from a sur- 
vey of the real situation. 

Given a few simple rules, the home micro 
user can conjure up alternative futures for a 
host of different events, but the reliability of 
the results depend on the accuracy of the data 
fed into the program or the precision of the 
rules by which the game is played. This fact 
can be demonstrated by a simulation of the 
outcome of football matches. 

Every Saturday afternoon during the soc- 
cer season, millions of pools punters eagerly 
await the football results and the chance of a 
bumper first dividend. Sadly, all but a few are 
disappointed, and must wait for another week 
and another chance. If you don’t want to wait 
for a whole week, then key in the first 
program and try your luck as many times as 
you like: 


10 POKE 23658,0 

20 DIM a(55): DIM r$(4,14) 

30 BORDER 0: PAPER 0: INK 7: CLS 
50 PRINT AT 0,7; INVERSE 1; 

“□□TREBLE CHANCED □” 

60 PRINT : PRINT 

80 PRINT “□□HOW MANY SELECTIONS 
FROM 550 □ □ □□ MATCHES DO YOU 
WISH TO MAKEO □□□□□□□□ 
□ □□□□(8-16)” 

90 INPUT n 

95 IF n < 1 OR n> 16 THEN CLS : GOTO 80 
100 PRINT : PRINT “OWHAT ARE THE 
NUMBER OF MATCHESD QTHAT YOU 
WISH TO SELECT (1-55)” 

120 FOR k = 1 TOn 
130 INPUT a(k) 

140 NEXT k 
150 CLS 

160 PRINT FLASH 1; PAPER 2;AT 
1 0,6;“ □ MATCHES IN PROGRESS^” 

170 FOR v = 1 TO 2000: NEXT v 
180 CLS 

190 PRINT AT 0,10; INVERSE 1; 

“□□RESULTSnD”" 

200 LET r$(1) = “HOME WIN”: LET 
r$(2) = “AWAY WIN” 

210 LET r$(3) = “GOAL — LESS DRAW”: LET 
r$(4) = “SCORE DRAW” 

220 FOR i = 1 TOn 
230 LET y = RND‘1 
240 IF y< = .5 THEN LET z$ = r$(1) 

250 IF y>. 5 AND y< =.75 THEN LET 
z$ = r$(2) 

260 IF y > .75 AND y < = .9 THEN LET 
z$ = r$(3) 

270 IF y>. 9 THEN LET Z$ = r$(4) 

280 PRINT “MATCH NUMBER”;TAB 
14;a(i);TAB 18;z$ 

290 NEXT i 
300 PRINT : PRINT 

310 INPUT “DO YOU WANT ANOTHER GO 

?□”;!$ 

320 IF t$ = “n” THEN GOTO 340 
330 GOTO 30 
340 STOP 

[E[E 

20 DIM A(55),R$(4) 


50 PRINT “□ M > || ^TREBLE CHANCE” 
80 PRINT“H||||HOW MANY SELECTIONS 
FORM 55 DO YOU”:PRINT“WISH TO 
MAKE” 

90 INPUT“(8 — 16)”;N 

100 PRINT “HUWHAT ARE THE NUMBERS 
OF THE MATCHES” 

110 PRINT “THAT YOU WISH 
TO”:PRINT“SELECT (1 —55)” 

120 FOR K = 1 TON 
130 INPUT A(K) 

140 NEXT K 

160 PRINT “□Jfl > U j3 MATCHES IN 
PROGRESS” 

170 FOR V = 1 TO 2000: NEXT V 
190 PRINT “□ a > II 3 RESULTS” 

200 R$(1) = “HOME WIN”: 

R$(2) = “AWAY WIN” 

210 R$(3) = “GOAL — LESS DRAW”: 

R$(4) = “SCORE DRAW” 

220 PRINT “HMATCH”:FOR 1 = 1 TO N 
230 Y = RND(1) 

240 IF Y < =0.5 THEN Z$ = R$(1) 

250 IF Y>0.5 AND Y< =0.75 THEN 
Z$ = R$(2) 

260 IF Y>0.75 AND Y< = 0.9 THEN 
Z$ = R$(3) 

270 IF Y> 0.9 THEN Z$ = R$(4) 

280 PRINT A(I)TAB(6)Z$ 

290 NEXT I 

310 PRINT" a DO YOU WANT ANOTHER GO 

(Y/N)?”; 

320 GET T$:IF T$ = “N” THEN 
PRINT“Q”:END 
330 IF T$ = “Y” THEN RUN 
340 GOTO 320 

B 

20 DIM A(55),R$(4) 

50 M0DE1 :PRINTTAB(1 4,3) 

“TREBLE CHANCE” 

80 PRINT”“HOW MANY SELECTIONS FROM 
55 MATCHES DODOYOU WISH TO 
MAKE (8-16)?” 

90 INPUT N 

100 PRINT“WHAT ARE THE NUMBERS OF 
THE MATCH ES”'“TH AT YOU WISH TO 
SELECT (1-55)?” 

120 FOR K = 1 TO N 
130 INPUTA(K) 

140 NEXT 
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■ 

RULES OF THE GAME 

■ 

GENERATING RANDOM NUMBERS 

■ 

SCORING 

■ 

TYPES OF RANDOM NUMBER 

■ 

DECIDING RESULTS 

■ 

SAMPLES AND SURVEYS 

■ 

SELECTING NUMBERS FROM 

■ 

SAMPLING 


A HAT OR BOX 

■ 

SINGLE-PASS SEARCH 


160 PRINTTAB(1 1 )“MATCHES 190 CLS:PRINTTAB(1 7,3) 

IN PROGRESS” “RESULTS” 

170 D = INKEY(500) 200 R$(1) = “HOME WIN”: 



R$(2) = “AWAY WIN” 

210 R$(3) = “GOAL — LESS DRAW”: 

R$(4) = “SCORE DRAW” 

220 FOR 1 = 1 TO N 
230 Y = RND(1) 

240 IF Y< = .5 THEN 
Z$ = R$(1) 

250 IF Y > .5 AND Y < = .75 THEN 
Z$ = R$(2) 

260 IF Y > .75 AND Y< = .9 THEN 
Z$ = R$(3) 

270 IF Y>. 9 THEN 
Z$ = R$(4) 

280 PRINT“MATCH NUMBER 
A(I)TAB(20)Z$ 

290 NEXT 

310 INPUT"“DO YOU WANT ANOTHER 
GO”,G$ 

320 IF G$ = “N” THEN END 
330 GOTO 50 


ES 


20 DIM A(55) 

30 CLS 

50 PRINT@9, “treble chance” 

60 PRINT:PRINT 

80 INPUT“HOW MANY SELECTIONS FROM 
55 □ □ □ □□ MATCHES DO YOU WISH 
TO MAKE □ □ □ □ □ (8 — 1 3)”;N 
90 IF N <8 OR N> 13 THEN 80 
100 PRINT“WHAT ARE THE NUMBERS OF 
THED □ □ □□ MATCHES THAT YOU 
WISH TO SELECT (1-55)11]” 

110 T = 0 

120 FORK = 1 TON 
130 INPUT A(K) 

140 NEXT K 
150 CLS 

160 PRINT@262, “matches in progress” 

170 FOR V=1 TO 2000: NEXT V 
180 CLS 

190 PRINT@1 2, “results” 

220 FOR 1 = 1 TO N 
230 Y = RND(0) 

240 IF Y< = .5 THEN Z$ = “HOME WIN” 

250 IF > .5 AND Y < = .75 THEN Z$ = “AWAY 
WIN” 

260 IF Y>. 75 ANDY <=.9 THEN 
Z$ = “GOAL -LESS DRAW” 
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270 IF Y > .9 THEN Z$ = 

“SCORE DRAW” 

280 PRINT“MATCH NUMBER 

A(I);TAB(17);Z$ 

290 NEXT I 
300 PRINT:PRINT 

310 INPUT“DO YOU WANT ANOTHER GO 

(Y/N)D”;T$ 

320 IF T$ = “N” THEN END 
330 GOTO 30 

There are no first dividends in this football 
simulation, but it may sharpen your feel for 
the real thing. In any event, it will 
demonstrate just how small is the chance of 
getting 24 points. 

RUN the program, and choose a number of 
matches between 8 and 16 in the range 1 to 55. 
This is equivalent to placing an X against the 
games chosen as likely draws. The results for 
the games you selected are then displayed on 
the screen. Scoring three points for a score 
draw, two for a no-score draw, one-and-a-half 
for an away win and one for a home win, the 
best eight matches in a typical selection would 
yield a total of only about 15 points — not a 
good start. 

In January and February, bad weather 
sometimes forces matches to be cancelled. 
When this happens, a panel of football experts 
decides how the results would have turned 
out if the games had actually taken place. In 
fact, both the program and the panel are 
simulators — they provide a symbolic repre- 
sentation of a real process. 

The program is not particularly difficult. 
After matches have been selected (Lines 120 
to 140) the main part of the simulation (Lines 
220 to 290) decides the outcome of each 
match, and prints the results. 


DECIDING RESULTS 


The basis on which the computer decides to 
print ‘score draw’, ‘no-score draw’, or another 
result, is the old pools punter’s rule — that a 
half of the matches played turn out to be home 
wins, a quarter away wins, and the remaining 
quarter result in draws. It has been assumed 
that the number of ties is divided between 
score draws and no-score draws in the ratio 
two to three. This ratio gives factors of f and f. 
When multiplied by J, the values are 0.1 and 
0.15, so that the probability (see pages 694 to 
700) of a score draw is 0.1. 

It is worth noting that previously, each 
match had only three possible results (home 
win, away win and draw). Then, the total 
number of different possible outcomes of 
eight matches were three to the power eight, 
or 6561. The introduction of a score draw 
category increased this dramatically to (4|8) 




or 65,536. If the chance of choosing a score 
draw is 1 in 10 (0.1), then the odds against 
picking eight score draws in eight matches is 1 
to 100 million, which is most discouraging. 

Choosing numbers in the treble chance is 
similar to deciding who wins a raffle. Usually, 
the tickets are placed in a large box or a hat, 
and someone, without looking, picks out the 
lucky number. The allocation of horses in 
the office Derby sweepstake can be organ- 
ized in a similar fashion, and the method 
is termed top hat simulation. 

Suppose a piece of paper is cut into 
four quarters. On the first quarter, the 
words ‘away win’ are written. ‘Draw’ 
is written on the second quarter and 
the other two pieces of paper bear 
the message ‘home win’. Now, if the 
four pieces of paper are folded and 
placed in a top hat, picking a piece of 
paper at random is equivalent to decid- 
ing the outcome of a football match. 

The process would be just as effective 
if one each of the numbers 1, 2, 3 and 
4 had been written on the four pieces of 
paper. The outcome could then have been 
decided by referring to a table: 


Number drawn Result 

1, 2 home win 

3 away win 

4 draw 


Similarly, when score draws are taken into 
account, if 40 pieces of paper bearing one each 
of the numbers 1, 2, ... 40 are placed in the 
top hat, the result of picking out a particular 
piece of paper can be interepted according to 
the following table: 


Number drawn 

1 to 20 
21 to 30 
31 to 36 
37 to 40 


Result 

home win 
away win 
no-score draw 
score draw 


The BASIC function RND, which generates 
random numbers, gives a computation 
equivalent to picking numbers out of a hat. 
You can specify that RND generates a number 
which is equally likely to take any value 
between 0 and 1. It only remains to rewrite 
the table in decimal format: 


Value of RND 

from 0 to 0.5 
greater than 0.5 to 0.75 
greater than 0.75 to 0.90 

greater than 0.90 to 1.00 


Result 

home win 
away win 
no-score 
draw 

score draw 


These are the actual values that are coded at 
Lines 230 to 270 in the program above. 


The treble chance simulation is only as 
good as the probability estimate used in the 
program. You could collect your own data 
from the Sunday papers, for example, and 
decide whether a greater proportion of score 
draws, say, is more usual. 


GENERATING RANDOM NUMBERS 


Many years ago, manual or mechanical 
methods, such as dice throwing, card- 
shuffling or the spinning of a roulette wheel, 
were used to generate random numbers. 
These methods were slow and tedious, so the 
famous American mathematician, John Von 
Neuman, proposed a mid-square technique. 
Starting with a four-digit number (the seed), 
the next ‘random number’ is obtained by 
multiplying the seed by itself, and then taking 
the middle four digits. For example, suppose 
the seed is 5272. Then the second number can 
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be generated by taking the middle four digits 
of (5272J2) or 27,793,984. The answer (7939) 
is, practically, random. A second number can 
be obtained by squaring 7939, and so the 
process continues. 

How, you may ask, can any mathematical 
process — which must be repeatable — produce 
true random numbers? The answer is that it 
can’t. The numbers so obtained, however, 
behave as though they are random, and are 
usually referred to as pseudo-random or 
quasi-random. Numbers generated by using 
RND are, in fact, pseudo-random. Unfortu- 
nately, the mid-square technique is not very 
useful for generating random numbers by 
computer. 

Besides being slow, the sequence quickly 
repeats, and once a zero is obtained, the whole 
process ends. Most home micros use a con- 
gruence method which uses remainders to 


generate pseudo-random sequences. The next 
program uses a simple formula and the I NT 
(integer) function to provide an example of 
this: 


20 BORDER 0: INK 7: PAPER 0: CLS 
30 PRINT AT 0,4; INVERSE 1; 

“□PSEUDO-RANDOM NUMBERS^”" 

40 INPUT “DHOW MANY NUMBERS 
?D”;n: LETs = 0 
50 LET x = .677829 # PEEK 23673/50 
60 LET x = x # 1 842.95 
70 LET x = x — INT (x) 

80 LET p = INT (x*1 000): PRINT 

“□□□□□”;.001’p, 

90 LET s = s + 1 
110 IF s< =n THEN GOTO 60 
120 STOP 

IE 

30 PRINTQ a > || fl PSEUDO 
-RND NUMBERSH” 

40 PRINT“H0W MANY NUMBERS”: 
INPUT N 
50 X = .677829'TI/60 
60 X = X*1 842.95 
70 X = X — INT (X) 

80 P = INT(X’1000):PRINT ,001'P, 

90 S = S + 1 

110 IF S< =N THEN 60 


30 M0DE1 :PRINTTAB(8,3) 

“PSEUDO -RANDOM NUMBERS” 

40 INPUT“H0W MANY NUMBERS”, N: 

S = 0 

50 X = ,677829*(TIME/1 00 — INT 
(TIME/100)) 

60 X=X*1 842.95 
70 X = X — INT(X) 

80 P = INT(X*1000):PRINT.001 *P; 

90 S = S + 1 

110 IF S< =N THEN 60 

sail 

20 CLS 

30 PRINT@3, “pseudo-random numbers” 

40 INPUT“H0W MANY NUMBERSD”; 

N:S = 0 

50 X = .677829*TIMER/50 
60 X = X f 1 842.95 
70 X = X — INT(X) 

80 P = INT(X*1000):PRINTLEFT$ 
(STR$(P/i000)+ <c nnnn 

□ □□”, 8 ); 

Q(7) C _ C . 1 

110 IF S< =N THEN 60 
120 END 

Line 50 uses the computer time function to 


set a different seed value for each run. The 
number .677829 is a fairly arbitrary constant. 
After the start value has been multiplied by 
another constant (Line 60), the remainder or 
decimal part is obtained. Change the value of 
the constants at Lines 50 and 60, and RUN the 
program again to see what son of results you 
obtain. 

For most purposes, it is wise — and far 
simpler — to use the RND function provided 
by your own computer. Varying the value of x 
in the expression RND(x) will usually enable 
you to select a repeatable sequence which is 
good for testing or to re-seed each time. 

Perhaps the most important point to re- 
member when writing programs is that the 
function RND provides a random variable and 
notan algebraic variable. RND(1 ) — or RND(0) 
on the Dragon and Tandy — listed at one line 
of a program will not take the same value as 
RND(1) listed elsewhere. 


SAMPLES AND SURVEYS 


Political opinion polls, consumer research 
organizations and governments use com- 
puters to generate random samples. It is 
important, for example, that a market re- 
search company interviewing 1000 people 
gets a typically varied sample, representative 
of a much larger group of the population. It 
would not do, for example, to choose people 
who were all members of a vintage car club, if 
the survey was supposed to check the national 
pattern of car ownership. 

The best way to determine that a sample is 
representative is to pick them randomly — this 
excludes any particular bias. 

This still holds true, even when you are 
working with a smaller sample who do have a 
particular common interest, like the members 
of a vintage car club, or the readers of 
INPUT. Selecting randomly is still import- 
ant here, if, say, you wanted to determine the 
spread of computer ownership of particular 
computers amongst INPUT readers. 

You may even wish to pick a sample from 
your own computerized club membership 
list. Basically, the sampling process is similar 
to a top-hat simulation in which several pieces 
of paper are picked from the hat. When 
simulating, however, each piece of paper is 
returned to the hat before the next is chosen. 
In sampling, once a piece of paper has been 
selected from the top hat, it is placed to one 
side before the next is chosen. Enter the next 
program to see how RND can be used to 
generate a random sample: 


10 DIM b$(10,16) 
20 DIM a$(10,16) 
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30 BORDER 0: PAPER 0: INK 7: CLS 
50 PRINT AT 0,7; INVERSE 1; 

“□RANDOM SAMPLINGD” 

90 PAUSE 100: CLS 

100 FOR i = 1 TO 10: READ a$(i): NEXT i 
110 INPUT “DSAMPLE SIZE ?D”;n 
120 FOR v=1 TO 10: LET b$(v) = a$(v): 

NEXT v 

130 FOR j = 1 TOn 

140 LET r=1 + INT (RND'1 0) 

150 IF b$(r) = “□□□□□□□□ 

□ □□□□□□□” THEN GOTO 140 
160 PRINT b$(r) 

170 LET b$(r) = 

180 NEXT j 

190 INPUT “□ ANOTHER SAMPLE (y/n) 

?D”;g$ 

200 IF g$ = “y” THEN CLS : GOTO 110 
210 STOP 

220 DATA“BONN”, “COPENHAGEN’ 1 , 
“LONDON” 

230 DATA “MADRID”, “MOSCOW”, “NEW 
YORK” 

240 DATA “PARIS”, “ROME”, “STOCKHOLM”, 
“VIENNA” 

BE [E 

20 DIM A$(10) 

50 PRINT “□H> U aRANDOM 
SAMPLINGS S” 

100 FOR 1 = 1 TO 10:READ A$(l): 

NEXT I 

110 INPUT “SAMPLE SIZE”;N:PRINT 
120 FOR V = 1 TO 10:B$(V) = A$(V): 

NEXT V 

130 FOR J = 1 TO N 
140 R = 1 + INT(RND(1)‘10) 

150 IF B$(R) = “” THEN 140 
160 PRINT B$(R) 

170 B$(R) = “” 

180 NEXT J 

190 PRINT “S3ANOTHER SAMPLE 

(Y/N)?” 

200 GET G$:IFG$ = “Y” THEN RUN 
210 IF G$< >“N” THEN 200 
215 PRINT“Q” 

220 DATA BONN, COPENHAGEN, LONDON 
230 DATA MADRID, MOSCOW, NEW YORK 
240 DATA PARIS, ROME, STOCKHOLM, 

VIENNA 


a 

20 DIM A$(10),B$(10) 

50 MODE1 :PRINTTAB(13,3)“RANDOM 
SAMPLING” 

90 D = INKEY(300) 

100 FOR 1 = 1 TO 10:READ A$(l): 
NEXT 

110 INPUT”“SAMPLE SIZE 1-1 0”,N 
120 FOR V = 1 TO 10:B$(V) = A$(V): 
NEXT 


130 FOR J = 1 TON 
140 R = RND(10) 

150 IF B$(R) =“” THEN 140 
160 PRINTB$(R) 

170 B$(R) = “” 

180 NEXT 

190 INPUT"“ANOTHER SAMPLE (Y/N)”,G$ 
200 IF G$ = “Y” THEN 110 
210 END 

220 DATA BONN, COPENHAGEN, LONDON 
230 DATA MADRID, MOSCOW, 

NEW YORK 

240 DATA PARIS, ROME, STOCKHOLM, 
VIENNA 




40 CLS 

50 PRINT@8, “random sampling”: 
PRINT:PRINT 

100 RESTORE:FOR 1 = 1 TO 10: 

READ A$(I):NEXT I 
110 INPUT“SAMPLE SIZED”;N: 
PRINT 

115 IF N>10 THEN 40 
120 FOR V = 1 TO 10:B$(V) = A$(V): 
NEXT V 

130 FOR J = 1 TON 



MAKING PREDICTIONS 


The simulation techniques described in 
this article have many of the elements you 
would need to devise a method of predict- 
ing results of a game or to carry out a 
survey. So it should be easy for you to 
modify the programs, or use certain sec- 
tions of them, in your own code. You 
could, for example, write a program to 
decide which matches will result in a score 
draw, then select eight of these matches at 
random and compare your predictions 
with actual results. This would require 
sections of the ‘treble chance’ and the 
‘random sampling’ programs, together 
with a few lines to link them and keep track 
of your selections. One use for such a 
program would be to remove the tendency 
for a person to choose matches on the 
coupon according to the separation of the 
numbers — they spread their choice at 
almost regular intervals on the coupon. In 
fact it could be interesting to compare both 
methods of choosing which of the games 
will produce the appropriate results. 


140 R = 1 +INT(RND(0) # 10) 

150 IF B$(R) = “” THEN 140 
160 PRINTB$(R) 

170 B$(R) = “” 

180 NEXT J 

190 PRINT:INPUT“ANOTHER SAMPLE 
(Y/N)D”;G$:PRINT:PRINT 
200 IF G$ = “Y” THEN 110 
210 END 

220 DATA BONN, COPENHAGEN, 

LONDON 

230 DATA MADRID, MOSCOW, 

NEW YORK 

240 DATA PARIS, ROME, STOCKHOLM, 

VIENNA 

After the data have been read in (Line 
100), a random integer R between 1 and 10 is 
generated (Line 140). Your computer will 
then print the Rth item (Line 160) on the list. 
Once an item has been selected, it must be 
removed from the data bank so that it is not 
selected a second time. Line 170 looks after 
that job. Notice that you could use RND(1 0) to 
generate a random variable in the range 1 to 
10, but the method shown at Line 140 follows 
on from the use of decimal fractions to test the 
range of probabilities. 

Usually, the data population from which 
the sample is drawn is far more extensive than 
the ten items in this example. For sampling 
such masses of data, the program above would 
be slow and inefficient. Suppose, for instance, 
a pollster wishes to pick 200 names from an 
electoral roll of 60,000 voters in a particular 
constituency, a program such as the one above 
would need to search the entire list of electors 
200 times. To avoid the long wait that this 
would entail, it would be better to use a 
single-pass method. 


THE SINGLE-PASS SEARCH 


The single-pass method would read through 
the electoral roll once from top to bottom. As 
each name is considered, a decision is made 
whether to include that person in the sample. 
This method can be programmed easily — 
make the following changes to the last 
program and reRUN it: 


50 PRINT AT 0,1; INVERSE 1; 

“□SINGLE -PASS RANDOM 
SAMPLINGD” 

120 LET a = n: LET c = 1 0 
130 FOR j = 1 TO 10 
140 IF a = 0 THEN GOTO 190 
150 IF RND'1 < =a/c THEN PRINT a$(j): 
GOTO 170 

160 LET c = c — 1 : GOTO 180 
170 LET a = a — 1 : LET c = c — 1 







KH 

50 PRINT “□ a > || fl SINGLE - 
PASS ” 

60 PRINT “U || a RANDOM 

SAMPLINGSBS” 

120 A = N:C = 10 
130 FOR J = 1 TO 10 
140 IF A = 0 THEN 190 
150 IF RND(1J< =A/C THEN PRINT 
A$(J):GOTO 170 
160 C = C — 1: GOTO 180 
170 A = A — 1:C = C — 1 

B 

50 MODE1 :PRINTTAB(5,3)“SINGLE PASS 
RANDOM SAMPLING” 

120 A = N:C = 10 
130 FOR J = 1 TO 10 
140 IF A = 0 THEN 190 


150 IF RND(1) < = A/C THEN 
PRINTA$(J):A = A — 1 
1 60 C = C — 1 

Also delete Line 170 


E 


50 PRINT@2, “single-pass random 
sampling”:PRINT:PRINT 
120 A = N:C = 10 
130 FOR J = 1 TO 10 
140 IF A = 0 THEN 190 
150 IF RND(0) < = A/C THEN 
PRINTA$(J):GOTO170 
1 60 C = C — 1 :GOTO 1 80 
170 A = A — 1:C = C — 1 

If now you enter 3, to select three items from 
the list, the first item (Bonn) is considered 
first. Should the function RND (Line 150) be 


less than 3/10, Bonn is selected. Copenhagen 
is considered next. If Bonn is already in the 
selection, Copenhagen will have a chance of 
being selected only if the function RND 
generates a value of less than 2/9. On the other 
hand, if Bonn is not already selected, 
Copenhagen’s chance will go up to 3/9. Lines 
160 and 170 update the probabilities. When 
you compare the results of selections from 
this program with those from the previous 
one, you should notice that the single-pass 
method gives samples in alphabetical order. 

At first sight, it might appear that with a 
list of only ten cities, the number of possible 
samples is small. This is not so. In fact, there 
are 120 possible different samples of size 
three, and 252 samples of size five. In a future 
article, you can see how thes'e ideas can be 
developed and used in a type of simulation 
referred to as modelling. 
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A surprising variety of patterns can 
be produced using very simple 
graphics routines that superimpose 
curves or build up patterns of dots. 
Here are a few ideas to try 


The way in which computers can be used to 
plot the orbit or trajectory of an object falling 
under gravity was described in the articles on 
pages 740 to 747 and 797 to 803. Such 
programs illustrate the simplest aspects of the 
old science of dynamics. Orbits can, however, 
generate patterns much richer and more 
interesting than the parabolas, circles and 
ellipses in which projectiles and planets move. 
This article will explain how some of these 
patterns can be produced by very elementary 
programming and graphics routines. 

Nowadays dynamics is once again a 
rapidly developing field of research. One 
reason for this is the realization that there are 
mathematical ideas underlying dynamics that 
can be applied much more widely than simply 
to the motion of bodies acted on by forces. 
Optical scientists interested in the deviation 
of starlight by refraction in the atmosphere, 
industrial chemists studying the progress of a 
reaction, and biologists concerned with the 
growth of populations of competing species, 
all find themselves using the mathematics of 
dynamics. Another reason is that 
computers — by enabling simple operations to 
be repeated many times— have led to the 
discovery of structural complexity often un- 
suspected on the basis of the restricted calcul- 
ations previously possible. A pattern or curve 
may need to be drawn many, many times 
before any structure begins to appear. The 
programs below show how this can be done. 

All the programs for the Commodore need 
a Simons’ Basic cartridge or INPUTs hi-res 
utility. The programs for the Vic need a 
Super Expander cartridge. 


FAMILIES OF ORBITS 


When assembled into a collection or family, 
orbits that are individually simple can display 
striking patterns. To see this for parabolas, 
enter and RUN the first program. 



5 LET A$ = FOR N = 1 TO 64: LET 
A$ = A$ + “D”: NEXT N 
10 BRIGHT 1: BORDER 0: INK 5: PAPER 0: 
CLS 

20 LET NMAX = 81 
30 LET DELT = .05 






■ 

USING PATTERNS IN SCIENCE 

■ 

PATTERNS OF DOTS 

■ 

FAMILIES OF ORBITS 

■ 

DOT CURTAINS AND THE RINGS 

■ 

ORBITS AND ENVELOPES 


OF SATURN 

■ 

FOCUSING LINES AT CUSPS 

■ 

COMPLETE CHAOS 

■ 

CATASTROPHE THEORY 

■ 

FISH POPULATIONS 



40 LET SX = 1 70/SQR 3: LET SY = 175 
50 FOR N = 1 TO NMAX 
60 LET A = PI*( — 1 + 2*N/NMAX) 

70 PLOT 128,80 

80 FOR T=0 TO 3 STEP DELT 

90 LET X = T'COS A: LET Y = T*(SIN A-T/2) 

100 IF Y< = — .4 THEN GOTO 120 

110 LET DX = SX'X + 1 28: LET 
DY = SY'Y + 80 

111 IF DX<0 OR DX > 255 OR DY<0 OR 
DY >175 THEN GOTO 120 

1 1 5 DRAW DX - PEEK 23677, DY - PEEK 23678 
117 PRINT AT 19,0;A$ 

120 NEXT T 
130 NEXT N 

BE 

10 HIRES 6,3:CO LOUR 6,3 
15 BLOCK 0,160,319,199,1 
20 NM = 81 
30 DE = .05 

40 SX = 1 60/SQR(3) :SY = 200 
50 FOR N = 1 TO NM STEP 2 
60 A = ti*( — 1+2‘N/NM) 

70 XX = 1 60: YY = 1 00 

80 FOR T = 0 TO 3 STEP DE 

90 X = T*COS(A):Y = T*(SIN(A) — T/2) 

100 IFY> — .4THENLINE XX,YY,SX*X + 1 60, 

1 00 - SY*Y,1 :XX = SX’X + 1 60:YY = 

100 — SY*Y:GOTO110 
105 T=3 
110 NEXT T,N 
130 GOTO 130 

KS 

10 GRAPHIC 2:COLOUR 6, 1,1,1 
15 DRAW 1,0, 000 TO 1023,800: 

PAINT 1,0,808 
20 NM = 81 
30 DE = .05 

40 SX= 51 2/SQR(3):SY = 1 000 
50 FOR N = 1 TO NM STEP 2 
60 A = 7i*( — 1 +2*N/NM) 

70 POINT 0,512,512 

80 FOR T = 0 TO 3 STEP DE 

90 X = T‘ COS(A) :Y = T* (SI N (A) - T/2) 

100 IFY> — ,4THEN:DRAW 1 TO SX'X + 512, 
51 2 — SY*Y:G0T01 1 0 
105 T=3 
110 NEXT T,N 
130 GOTO 130 


B 

10 MODE0 
20 NMAX = 81 
30 DELT = .05 

40 SX = 800/SQR(3):SY = 1 300 
50 FOR N = 1 TO NMAX 
60 A = Pl‘( — 1 + 2*N/NMAX) 

70 MOVE 640,341 
80 T=0: REPEAT 

90X = rCOS(A) : Y = T* (SIN (A) - T/2) 

100 IF Y> -0.4 THEN DRAW SX*X + 640, 
SY*Y + 341 

110 T=T + DELT:UNTIL Y< = — .4 ORT>3 
120 NEXT 


ES 


10 PMODE4,1:PCLS1:SCREEN1,0 
20 NM = 81 :PI = 4*ATN(1 ) 

30 DE = .05 

40 SX = 1 60/SQR(3):SY = 230 
50 FOR N = 1 TO NM 
60 A = Pl*( — 1 + 2 # N/NM) 

70 DRAW“BM127,118” 

80 FOR T = 0 TO 3 STEP DE 

90 X = T*COS(A):Y = T*(SIN(A) — T/2) 

100 IF Y> -.4 THEN LINE — (SX*X -h 1 27, 

1 18 -SY # Y), PRESET ELSE T = 3 
110 NEXT T 
120 NEXT N 
130 GOTO 130 

This program simulates the paths of falling 
drops of water sprayed from the head of a 
garden sprinkler, or, in a more modern 
application, neutrons squirted from a thin 
pipe connected to a reactor. In this case the 
family of orbits consists of all the parabolic 
paths that emerge from a particular point in 
different directions but with the same speed. 
The pattern formed by this family is the outer 
curve which each orbit touches. This curve, 
called the envelope of the family, happens, 
itself, to be a parabola in this case (in gunnery 
it is called the ‘bounding parabola’ because it 
is the boundary of the region that can be 
reached by projectiles of a fixed initial speed). 
It is important to realize that the envelope is a 
property of the whole family of parabolic 
orbits and has no meaning for any single one 
of them. Thus envelopes are perfect ill- 
ustrations of the fact that the whole can be 
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greater than the sum of its parts. 

In the program, the equation for the orbits 
appears on Line 90. X is horizontal distance, Y 
is vertical distance, T is time (starting from 
zero at the instant of emission), and each orbit 
in the family is labelled by A, which is an angle 
giving the direction in which it starts out. 
There are NMAX or NM such directions; try 
changing the value of NMAX or NM in Line 20. 


FOCUSING 


Even straight lines can form families with 
interesting envelopes, as the second program 
shows. 



10 BORDER 0: INK 7: PAPER 0: CLS: 

LET N = 2 

30 LET Y0 = 80/N/N 

40 PLOT 0, — 6— Y0‘2 

50 FOR X = 0 TO 255 STEP 2 

60 LET Y = 169 — Y0 — Y0'COS (N*2*PI*X/255) 


70 DRAW X - PEEK 23677,1 75 - Y - PEEK 
23678 

80 LET XT = X + 2'Y0'N‘PrY/255’SIN (N‘2 1 
PI’X/255) 

85 LET Y1 =0: IF XT <0 THEN LET XT = 0: 
LET Y1 =Y + 255'X/(2'N'PI‘Y0‘SIN (N* 
2'PI'X/255)) 

86 IF XT > 255 THEN LET XT = 255: LET 
Y1 =Y — (255 — X)'255/(2*N'PrY0'SIN 
(N'2'PrX/255)) 

90 DRAW XT - PEEK 23677,1 75 -Y1- PEEK 
23678 

100 PLOT X,1 75 — Y 

110 NEXT X 

120 GOTO 120 

KH 

10 HIRES 6,3:COLOUR 6,3 

20 N = 2 

30 Y0 = 80/N/N 

40 XX = 0:YY = INT(1 60 — Y0*2) 

50 FOR X = 0 TO 319 STEP 2 


60 Y = 1 60 - Y0 - Y0'COS(N‘2*7t'X/31 9) 
70 LINE XX,YY,X,Y,1 :XX = X:YY = Y 
80 XT = X + 2*Y0 *N*tt*Y/31 9*$IN(N*2V 
X/319) 

85 Y1 =0:IF XT<0 THEN XT = 0:Y1 =Y + 
1 99*X/(2*N*7r*Y0*SIN(N*2*7t*X/1 99)) 

86 IF XT >31 9 THEN XT=319:Y1 = Y- 
(199 — X)*199/(2'N'7t‘Y0'SIN(N*2*7t’ 
X/199)) 

90 LINE XX,YY,XT,Y1 ,1 
110 NEXT X 
120 GOTO 120 

IE 

10 GRAPHIC 2:C0L0R 6, 1,1,1 
20 N = 2 

30 Y0 = 200/N/N 
40 XX=0:YY = INT(512 — Y0‘2) 

50 FOR X = 0 TO 1023 STEP 20 
60 Y = 51 2 — Y0 — Y0*COS(N*2*7tr*X/1 023) 
70 DRAW 1,XX,YYT0 X,Y:XX = X:YY = Y 
90 XT = X + 2*Y0*N VY/1 023*SIN(N*2 *tt* 
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X/1023) 

85 Y1 =0:IF XT<0 THEN XT = 0:Y1 =Y + 
1023'X/(2'N'7r'Y0'SIN(N'2VX/1023)) 

86 IF XT >1023 THEN XT = 1023:Y1 = Y- 
(1023 — X)'1023/(2'N'7r'Y0'SIN(N'2V 
X/1023)) 

90 DRAW 1,XX,YYT0 XT,Y1 
110 NEXT X 
120 GOTO 120 

a 

10 MODE0 
20 N = 2 

30 Y0 = 400/N/N 

40 MOVE0,Y0*2 + 30 

50 FOR X = 0 TO 1283 STEP 9 

60 Y = Y0 + 30 + Y0'COS(N'2'PI'X/1279) 

70 DRAW X,Y 

80 XTOP = X + 2'Y0'N'PI'(1023 - Y)/1 279' 
SIN(N'2’PI’X/1279) 

90 DRAW XTOP, 1023 
100 MOVEX,Y 
110 NEXT 


E 


10 PMODE4,1 :PCLS1 :SCREEN1 ,0 
20 N = 2: PI = 4*ATN(1 ) 

30 Y0 = 80/N/N 

40 DRAW"BM0,” + STR$(INT(1 86 - Y0'2)) 

50 FOR X = 0 TO 255 STEP 2 

60 Y = 1 86 — Y0 — Y0’COS(N'2'PrX/255) 

70 LINE — (X,Y),PSET 

80 XT=X + 2'Y0'N'PI'Y/255‘SIN(N'2'PI'X/ 
255) 

85 Y1 =0:IF XT <0 THEN XT = 0:YZ = Y + 
255*X/(2*N*PI*Y0‘SIN(N*2*Pt*X/255)) 

86 IF XT > 255 THEN XT = 255:Y1 = Y- 
(255 — X)*255/(2'N*PI*Y0*SIN(N*2*PI*X/ 
255)) 

90 LINE- (XT, Y1), PRESET 
1 00 DRAW“BM” + STR$(X) + + STR$ 

(INT(Y)) 

110 NEXT 
120 GOTO120 

This program simulates light rays bent by 
refraction through a wavy, curved surface, for 
example sunlight refracted by ripples on 
water in a swimming pool. The family of 
orbits consists of all the straight lines at right 
angles to a wavy curve of sine form. For light 
rays the envelope is the curve corresponding 
to focusing. Focusing is particularly intense 
near the troughs (minima) of the sine wave, 
where the envelope has sharp points called 
cusps. These cusps are the bright points of 
light you see amongst the ripples. The 
existence of cusps is predicted by the 
recently developed mathematics of envel- 
opes, called catastrophe theory (whose drama- 
tic name originated in applications of the 


same mathematics to the collapse of bridges 
and the capsizing of ships). 

In the program, the wavy initial curve is 
specified in Line 60 and the rays that start out 
from it are defined in Line 80. The number of 
troughs of the wavy curve is N; try changing 
the value of N in Line 20 (N = 1 is especially 
recommended). 

Trajectories that are themselves wavy 
sine curves are assembled into a family in 
the third program. 


10 BORDER 0: PAPER 0: INK 7: CLS 
20 LET QM = SQR (2'LN (3)) 

30 FOR Q= -QM TO QM'1.001 STEP QM/30 
40 PLOT 0,75 + Q'75/QM 
50 FOR T = 0 TO 3‘PI STEP .2 
60 LET X = Q'COS (EXP ( - Q'Q/2)'T) 

70 DRAW (T240/3/PI)- PEEK 

23677,75 + (X'75/QM) - PEEK 23678 
80 NEXT T 
90 NEXT Q 

SS 

10 HIRES 0,1 :C0 LOUR 1,6:MULTI 4,3,7 
20 QM = SQR(2'LOG(3)) 

30 FOR Q = - QM TO QM'1.001 STEP 
QM/30 

40 XX = 0:YY = INT(1 00 - Q'90/QM) 

50 FOR T=0 TO 3'n STEP .2 
60 X = Q‘C0S(EXP( — Q'Q/2)*T) 

70 LINE XX,YY,T'17/3'7r,100-X'90/QM,RND 
(1)-3 + 1 

75 XX = T*1 7/3*tt:YY = 100- X'90/QM 
80 NEXT T,Q 
90 GOTO 90 

EB 

10 GRAPHIC LCOLOR 6,1 ,3,7 
20 QM = SQR(2'LOG(3)) 

30 FOR Q = - QM TO QM'1.001 STEP 
QM/30 

40 POINT 0,0,INT(512 — Q‘500/QM) 

50 FOR T = 0 TO 3'n STEP .6 
60 X = Q'C0S(EXP( — Q*Q/2)'T) 

70 DRAW RND(1)'3 + 1 TO T*116/ 3*tt, 

512- X- 500/ QM 
80 NEXT T,Q 
90 GOTO 90 


□ 

10MODE0 

20 QM = SQR(2'LN(3)) 

30 FOR Q = - QM TO QM'1.001 STEP QM/30 
40 MOVE 0,Q'450/QM + 500 
50 FOR T = 0 TO 3* PI STEP .2 
60 X = Q‘C0S(EXP( — Q'Q/2)'T) 

70 DRAW T'l 200/3/PI, X'450/QM + 500 
80 NEXT 
90 NEXT 


ES 


10PMODE4,1:PCLS1:SCREEN1,0 
20 QM = SQR(2*LOG(3)):PI = 4*ATN(1 ) 

30 FOR Q= - QM TO QM'1.001 STEP 
QM/30 

40 DRAW“BM0,” + STR$(INT(1 00 - 0*90/ 
QM)) 

50 FOR T = 0 TO 3*PI STEP .2 
60 X = Q*COS(EXP( — Q*Q/2)*T) 

70 LINE - (F240/3/PU 00 - X f 90/QM), 
PRESET 
80 NEXT T 
90 NEXT Q 
100 GOTO 100 

This program simulates light rays passing 
along a thin glass optical fibre whose refrac- 
tive index varies across its width, or (on a 
much smaller scale) electrons winding be- 
tween planes of atoms in a crystal placed in 
the beam of an electron microscope. The 
family consists of orbits starting out parallel 
to each other at the left-hand edge of the 
screen; each orbit undulates about the central 
line at a rate that depends on its starting point. 
Although this family is very different from 
the last program, the envelope curves also 
display intense focusing at cusp catastrophe 
points. 

In the program, orbits are specified in Line 
60 by giving their distance X from the centre 
line at time T. Different orbits in the family 
are labelled by Q. There are 30 orbits in the 
program as written; to get more or fewer, 
change the number at the end of Line 30. 


DOT PATTERNS 


Assembling orbits into families is not the only 
method of obtaining interesting patterns. 
Another way is to follow the orbits for very 
long times , with the result that considerable 
complexity can develop even when the math- 
ematics of the orbit is relatively simple. In 
displaying these orbits it is not usually advis- 
able to plot the continuous curve giving the 
position at every instant, because over long 
times this would just fill the screen with a 
mess like tangled wool. Instead the orbit is 
plotted at regular intervals (once every sec- 
ond, for example), as though the trajectory 
were viewed in the flashing light of a strobos- 
cope. Examples of the resulting dot patterns 
are given in the next three programs. The 
position of dots on the screen does not always 
correspond to the real, spatial position of 
physical objects whose motion the programs 
simulate; sometimes it is an abstract repres- 
entation, in which horizontal screen position 
corresponds to position and vertical screen 
position to speed. 









The fourth program will take several 
minutes to RUN. 


10 BORDER 0: PAPER 0: INK 7: CLS 
30 LET A = 76.11 

40 LET ALF = A'PI/180: LET C = COS (ALF) 
50 LET S = SIN (ALF) 

60 LET N MAX = 200 
70 LET M = 52 
80 FOR J = 1 TO M 
90 LET X = 0: LET Y = J/M 
100 FOR N = 1 TO NMAX 
1 1 0 LET W = X 

120 LET X = X*C — (Y — X'X)"S: LET 
Y = WS + (Y-W*W)X 
1 30 IF ABS (X) > 4 OR ABS (Y) > 4 THEN 
GOTO 860 

135 IF X>1 OR Y > 1 THEN GOTO 150 
140 PLOT X*128 + 128,Y'85 + 85 
150 NEXT N 
160 NEXT J 

BB 

10 HIRES 0,1 :COLOUR 1,6:MULTI 5,3,7 
20 A = 76.1 

30 AL = A*7t/1 80:C = COS(AL) 

40 S = SIN(AL) 

50 NM = 200 

60 M = 52 

70 FOR J = 1 TO M 

80 X = 0:Y = J/M:CL = RND(1)‘3 + 1 

90 FOR N = 1 TO NM 

100 W = X 

1 1 0 X = X*C - (Y - X*X)*S:Y = W’S + (Y - 
W*W)X 

1 20 IF ABS(X) > 4 OR ABS(Y) > 4 THEN 1 60 
125 IF ABS(Y) > 1 OR ABS(X) > 1 THEN 140 
1 30 PLOT 80 + X'79,1 00 - Y'99,CL 
140 NEXT N,J 
150 GOTO 150 


IE 

10 GRAPHIC 1:COLOR 6, 1,3, 7 
20 A = 76.1 

30 AL = A'ti/1 80:C = COS(AL) 

40 S = SIN(AL) 

50 NM = 200 

60 M = 52 

70 FOR J = 1 TOM 

80 X = 0:Y = J/M:CL=RND(1)"3 + 1 

90 FOR N = 1 TO NM 

100W = X 

1 1 0 X = X*C - (Y - X’X)*S:Y = WS + 
(Y-WW)X 

1 20 IF ABS(X) > 4 OR ABS(Y) > 4 THEN 1 60 
125 IF ABS(Y) >1 OR ABS(X)>1 THEN 140 
1 30 POINT CL, 51 2 + X*51 2,51 2 - Y’51 2 
140 NEXT N,J 
150 GOTO 150 


a 

10 MODE0 
20 A = 76.11 

30 ALF = A*PI/180:C = COS(ALF) 

40 S = SIN(ALF) 

50 NMAX = 200 
60 M = 52 
70 FOR J = 1 TO M 
80 X = 0:Y = J/M 
90 FOR N = 1 TO NMAX 
100 W = X 

1 1 0 X = XX - (Y - X"X)’S:Y = WS + 
(Y-WW)X 

120 IF ABS(X) >4 OR ABS(Y) >4 THEN END 
1 30 PLOT 69,X‘640 + 640,Y*51 2 + 51 2 
140 NEXT 
150 NEXT 


LTH 


ID 


10PMODE4,1:PCLS1:SCREEN1,0 
20 A = 76.1 1 :PI = 4*ATN(1 ) 


30 AL = A # PI/1 80:C = COS(AL) 

40 S = SIN(AL) 

50 NM = 200 
60 M = 52 
70 FOR J = 1T0M 
80 X = 0:Y = J/M 
90 FORN = 1TONM 
100W = X 

1 1 0 X = X # C - (Y - X*X)*S:Y = W # S + 

(Y — W # W) # C 

1 20 IF ABS(X) > 4 OR ABS(Y) > 4 THEN 1 60 
125 IF ABS(Y) > 1 OR ABS(X) > 1 THEN140 
130 PRESET(1 28 + X # 1 28,96 — Y*96) 

140 NEXT N 
150 NEXT J 
160 GOTO160 

This program simulates motion of subatomic 
particles (such as protons or electrons) in an 
accelerator, or the windings of a line of force 
in the magnetic bottle of a fusion power 
device. Mathematically, what the program 
does is to take a series of initial points and 
move them about the screen by repeating a 
rule that says: ‘rotate about the middle of the 
screen by a fixed angle, apart from a slight 
modification’. Without the ‘slight modific- 
ation’ the orbits of the points would all be 
circles, and indeed those near the centre are 
roughly circles. But the modification has a 
dramatic effect on points initially far from the 
centre: their orbits may be a series of ‘islands’, 
or they may escape from the screen altogether. 
More refined computer graphics, using high 
magnification, reveals tiny islands every- 
where, distributed amongst the large ones. 

In the program, the number of initial 
points is M, specified in Line 60; if you get 
tired of waiting for the picture you can reduce 
this number. Each of these points is plotted 
for NMAX or NM repetitions of the transform- 
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ation rule; NMAX or NM is specified in Line 
50. The rule itself is contained in Lines 100 
and 110, which describe how the horizontal 
and vertical screen coordinates X and Y are 
altered at each repetition. The angle of the 
unmodified rotation, in degrees, is A, speci- 
fied in Line 20; you should experiment with 
different values of A (try 90). 

An unexpected pattern made of dots gen- 
erated by a single initial point is produced by 
the next program which you should now enter 
and run: 


20 LET K = 51 .3: BORDER 0: INK 7: PAPER 0: 
CLS 

30 LET X = 1/PI: LET P = 0 
40 LET A = 1/SQR 5 
50 FOR N = 1 TO 10000 


60 LET Y = X — .5: LET X = X + A — INT (X-f A) 
70 LET P = P — Y 
80 PLOT X’255,P # K + 85 
90 NEXT N 


BE 

10 HIRES 0,1:COLOUR 1,6: 

MULTI 5,3,7:K=60 
20X = 1/tt:P = 0 
30 A = 1/SQR(5) 

40 FOR N = 1 TO 10000 

50 Y = X — ,5:X = X + A — INT(X + A) 

60 P = P-Y 

70 PLOT X'159,100 — P‘K,RND 
(1)*3 + 1 
80 NEXT N 
90 GOTO 90 

BE 

10 GRAPHIC 1:C0L0R 6,1 ,3,7: 

K = 300 

20 X = 1/tc:P = 0 


30 A = 1/SQR(5) 

40 FOR N = 1 TO 10000 

50 Y = X — .5:X = X + A — INT(X + A) 

60 P = P — Y 

70 POINT RND(1)‘3 + 1,X'1023, 
512— P'K 
80 NEXT N 
90 GOTO 90 

B 

10 MODE0:K = 300 
20 X = 1/PI:P = 0 
30 A = 1/SQR(5) 

40 FOR N = 1 TO 10000 

50 Y = X — .5:X = X + A — INT(X + A) 

60 P = P — Y 

70 PLOT 69,X*1279,P*K + 512 
80 NEXT 




10 PMODE4,1:PCLS1:SCREEN1,0: 

K = 60 

20 X = 1/(4*ATN(1)):P = 0 
30 A = 1/SQR(5) 

40 FOR N = 1 TO 10000 

50 Y = X — .5:X = X + A — INT(X + A) 

60 P = P — Y 

70 PRESET (X*255,1 28 — P*K) 

80 NEXT N 
90 GOTO90 

This program is an abstract simulation of 
resonance , where the frequencies of two phys- 
ical effects may clash. For example, one 
frequency might be that of an asteroid’s 
motion round the sun, and the other might be 
the frequency with which the asteroid is 
disturbed by the gravitational pull of the 
planet Jupiter. In this case an important 
question is: do the disturbing pulls mount up 
and throw the asteroid out of its orbit, or do 
they force it into a stable orbit? The answer is 


that this depends on the ratio of the two 
frequencies (that is, on one divided by the 
other). The particles in Saturn’s rings are 
affected by similar forces. 

In the program, the ratio is called A and its 
value (which must be less than 1) is specified 
in Line 30. The number of repetitions of the 
resonance transformation is 10000 and is 
specified at the end of Line 40; reduce this 
value if you get tired of waiting. The trans- 
formation itself is on Lines 50 and 60, and is 
the rule for changing the horizontal (X) and 
vertical (P) positions of the dot on the screen. 

As written, the delicate curtain pattern is 
produced by A = 1/SQR(5) which equals 
0.4472136. This is not the ratio of two whole 
numbers (it is an ‘irrational number’) and 
corresponds to nonresonance. To get reson- 
ance you should try a number that is the ratio 
of two whole numbers, such as A = 9/20 
which equals 0.45. Then try a different 
irrational number, such as A = 1/PI. (If you 
want to reduce the vertical scale to display 
more of the picture, diminish the value of K. 


CHAOS 


It is a remarkable fact that some simple rules 
produce no pattern at all. Such an example is 
the program below which you should enter 
and RUN: 


10 BORDER 0: PAPER 0: INK 7: CLS 
20 LET X = 1/PI 
30 LET Y = 1/PI 
40 FOR M = 1 TO 10000 
50 LET X = X + Y — INT (X + Y) 

60 LET Y = X + Y — INT (X + Y) 

70 PLOT X'255,Y*175 
80 NEXT M 







BE 

10 HIRES 0,1:COLOUR 6,6: 

MULTI 5,3,7 
20 X = 1/7t 
30 Y = 1 /tt 

40 FOR M = 1 TO 10000 
50 X = X + Y — INT(X + Y) 

60 Y = X + Y — INT(X + Y) 

70 PLOT X'159,199— Y*199, 
RND(1)*3 + 1 
90 NEXT M 
90 GOTO 90 

KB 

10 GRAPHIC 1:COLOR 6, 1,3, 7 
20 X = 1/7t 
30 Y = 1/7r 

40 FOR M = 1 TO 10000 
50 X = X + Y — INT(X + Y) 

60 Y = X + Y— INT(X + Y) 

70 POINT RND(1)’3 + 1,X*1023, 
1023— Y*1 023 
80 NEXT M 
90 GOTO 90 

B 

10MODE0 
20 X = 1 /PI 
30 Y = 1/PI 

40 FOR M = 1 TO 10000 
50 X = X + Y— INT(X+Y) 

60 Y = X + Y — INT(X + Y) 

70 PLOT69,X*1 279X1023 
80 NEXT 


ES 


10PMODE4,1:PCLS1:SCREEN1,0 
20 PI = 4*ATN(1 ):X = 1 /PI 
30 Y = 1 /PI 

40 FORM = 1TO10000 
50 X = X + Y — INT(X + Y) 

60 Y = X + Y — INT(X + Y) 

70 PRESET (X*255,1 92 — Y*1 91 ) 

80 NEXT M 
90 GOTO 90 

This program simulates the erratic bouncing 
of metal spheres in a pinball machine, or the 
motion of molecules in a gas, or of a roulette 
wheel, or indeed any dynamic system whose 
orbits are so unpredictable as to be indistin- 
guishable from what would be generated by 
purely random processes. And yet the ten 
thousand points on the screen are not the 
result of random sprinkling but are generated 
by the program from a single point by 
repeating a rule that is purely deterministic — 
that is, it contains no random element. 

The rule is based on regarding the screen as 
the ‘unit square’ on which horizontal distance 


X and vertical distance Y range from zero to 
one, and operates in three stages, constituting 
Lines 50 and 60 of the program. First, the X 
and Y coordinates of a point are added to 
produce a new X; second, this new X is added 
to Y to produce a new Y; third, if either of the 
new X or Y lies outside the range 0 to 1, an 
appropriate whole number is subtracted in 
order to bring the transformed point back into 
the unit square (in the program this subtrac- 
tion is implemented by the I NT function in 
Lines 50 and 60). 


FISH POPULATIONS 


How does the population of fish in a pond 
change over many generations? This depends 
on the rule that determines how the popul- 
ation changes from one generation to the next. 
Such a rule must incorporate both the ten- 
dency of the population to increase because 
each set of parents produces more than two 
offspring, and the tendency of the population 
to decrease when it gets so large that the finite 
food supply in the pond cannot sustain it. 
Depending on the precise balance between 
these two dependencies a fish population may 
settle down to a stable fixed value, or alternate 
regularly between two or more values, or 
change apparently randomly between suc- 
cessive generations. 

The final program employs graphics and 
sound to illustrate these different possi- 
bilities. Horizontal screen position, denoted 
by A, corresponds to the balance between 
breeding and food supply and hence to the 
rule relating successive generations. Vertical 
screen position, denoted by Y, corresponds to 
the population, represented by a point jump- 
ing up or down at each generation. The 
population Y is plotted on the screen only 
when it has settled down to its stable value or 
set of alternating values, called the attractor 
set. But the way in which the population 
homes in on the attractor can be heard 
because the program makes the computer 
emit a sound whose pitch (frequency) is 
proportional to the population. 



10 BORDER 0: INK 7: PAPER 0: CLS 
20 LETS = .03 LET NMIN = 50: LET 
N MAX = 80 

30 FOR A = 2.8 TO 4 STEP S 
40 LET Y = 1/PI 
50 FOR N = 1 TO NMAX 
60 LET Y = A*Y*(1 — Y) 

70 IF N> NMIN THEN PLOT 255* 

(A — 2.8)/1.2,Y # 175 
80 BEEP .0075, Y # 20 
90 NEXT N 
100 NEXT A 


KB 

10 HIRES 0,1 :C0L0UR 1,6: 

MULTI 5,3,7 

15 POKE 54296,1 5:P0KE 54277,64 
20 S = 1/1 60:NN = 50:NX = 80 
30 FOR A = 2.8 TO 4 STEP S 
40 Y=.25/ATN(1) 

45 POKE 54276,33 
50 FOR N = 1 TO NX 
60 Y=A‘Y*(1 — Y) 

70 IF N>NN THEN PLOT 1 59*{A — 2.8)/1 .2, 
199 — Y*199,RND(1)*3 + 1 
80 POKE 54273,1 + 255'Y 
90 NEXT N 
95 POKE 54276,32 
100 NEXT A 
110 GOTO 110 

KB 

10 GRAPHIC 1:C0L0R 6,1 ,3,7 
15 POKE 36878,15 
20 S = .01 :NN = 50:NX = 80 
30 FOR A = 2.8 TO 4 STEPS 
40 Y = .25/ATN(1) 

50 FOR N = 1 TO NX 
60 Y = A*Y'(1 — Y) 

70 IF N>NN THEN:POINT RND(1)*3 + 1, 

1 023* (A — 2.8)/1 .2,1 023 — Y*1 023 
80 POKE 36876,1 28 + 127*Y 
90 NEXT N 
95 POKE 36876,0 
100 NEXT A 
110 GOTO 110 


B 

10MODE0 

20 S = .03:NMIN = 50:NMAX = 80 
30 FOR A = 2.8 TO 4 STEP S 
40 Y= 1/PI 

50 FOR N = 1 TO NMAX 
60 Y = A’Y’(1 -Y) 

70 IF N> NMIN THEN PLOT69,1279* 
(A — 2.8)/1 .2,Y*1 023 
80 S0UND1 , — 1 5,255*Y,1 
90 NEXT 
100 NEXT 


l t _*< 


Q 


10 PMODE4,1:PCLS1:SCREEN1,0 
20 S = 1/127:NN = 50:NX = 80 
30 FORA = 2.8 T04 STEPS 
40 Y = .25/ATN(1) 

50 F0RN = 1T0NX 
60 Y = A'Y'(1 — Y) 

70 IF N>NN THENPRESET(255'(A — 2.8)/1.2, 
192— Y‘191) 

80 S0UND1 +255'Y,1 
90 NEXT N 
100 NEXT A 
110 GOTO110 








It is clear that the attractor undergoes drastic 
changes as the amount of food is slowly 
altered. To start with (on the left of the 
screen), the attractor is a single point, indicat- 
ing an eventually stable fish population. This 
means there is enough food for the fish to live 
healthily and breed well. If the population 
increases too much there will be less food to 
go around so some fish will die of starvation. 
This leaves more food for the survivors so 
they grow and breed well again. Eventually 
the numbers in each generation settle down 
to a constant stable population. In the 
program you’ll hear the sound oscillating 
between each generation but settling down to 
a stable value each time. 

Suddenly, at a certain value of A, as the 


food is increased past a certain level, the 
attractor is seen and heard to branch into two 
or bifurcate , indicating a population eventu- 
ally alternating between two values. Later, 
with even more food, the attractor bifurcates 
again, indicating four alternating population 
values. More and more divisions occur, form- 
ing an infinite sequence of which only the first 
few are resolved by this program. The suc- 
cessive bifurcations accumulate at a finite A 
value, corresponding to a fish population that 
fluctuates among infinitely many values with- 
out ever settling down. 

Much excitement has been generated by 
the discovery that this ‘bifurcation tree’, 
ending in chaos, arises in a wide range of 
mathematical contexts. An important applic- 


ation is to the way in which the motion of a 
flowing liquid changes in stages from smooth 
(stable attractor) to turbulent (chaotic at- 
tractor) as the speed increases (from the lazy 
river to the raging torrent). The same thing 
can be seen as smoke from a cigarette rises 
gently then suddenly spirals and swirls into 
turbulent eddies. 

In this program the population evolution 
rule is specified in Line 60. Use of the sound 
command gives a vivid impression of how the 
orbit reaches the attractor, but does slow the 
program down. To get a clearer picture of the 
attractor itself in a reasonable time first delete 
Line 80. You can increase the resolution by 
changing Line 20 so that S = 0.005, NMIN or 
NN = 200 and NMAX or NX = 300. 
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TAILORING 

SPREADSHEETS 


Find out where your money goes or 
plan out the future of your business 
with this handy spreadsheet 
program. Add some more lines to 
the program started last time 


Since a spreadsheet starts out as a blank sheet 
of paper (or a blank screen!) it is sometimes 
difficult to know exactly what to use it for. 
The examples given last time and the ideas 
below should help, and you should be able to 
tailor one of them to your own needs. And of 
course, there is no need to limit yourself to 
just one sheet. The program can be used to set 
up any number of spreadsheets which can be 
saved and reloaded at any time. 

Here are a few examples. You could have 
one sheet to record and plan out your house 
expenses, where entries would go under head- 
ings such as rates, mortgage, insurance, re- 
pairs, etc. But if you were making a lot of 
repairs you might want to set up a separate 
sheet for these. You could itemise the differ- 
ent types of repairs, as well as decorations and 
furnishings for individual rooms, showing 
how much was spent each month or each 
quarter. The program could then total the 
expenditure for each type of repair over the 
whole house, or the total expenditure for each 
room. 

Yet another sheet could include all general 
family expenses such as food and drink, 
clothing, car, holidays and entertainment. It 
could list these under different weeks or 
months, or under different members of the 
family. 

But because the sheet starts off completely 
blank you can use it for any information that 
needs to be laid out in a logical way. For 
example a club membership spreadsheet 
could list members’ names, numbers and 
subscriptions paid, as well as listing and 
adding up the attendances each week— using 
labels for most of the entries. 

A spreadsheet is also ideal for entering 
details of any survey. In fact if you look back 
at the article on pages 269 to 275 you’ll see 
that the spreadsheet is really just a much more 
sophisticated version of a two-dimensional 
array. The spreadsheet allows you much more 
control over the way the entries are made and, 
because you can enter equations as well, the 
result of the survey can be calculated 
immediately. 

Business uses are more obvious, and are 
virtually unlimited. Spreadsheets can hold 
details of invoices, displaying details of the 
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items along with their cost, VAT, discounts 
and so on. They can work out staff salaries, 
where the names of the staff, hours worked, 
rate of pay, allowances and tax are all listed or 
worked out by the spreadsheet. They can be 
used for stocktaking, for general accounts and 
for budgeting — including making alternative 
forecasts for the future at the push of a button. 


ENTERING THE PROGRAM 


The section of program given here joins on to 
the one given last time. The remainder ap- 
pears in the next article which also gives 
detailed instructions on how to use it. So 
LOAD in the last part, add these lines and SAVE 
it ready for next time when you will be able to 
RUN the complete program. 

420 FOR a = fc TO tc: FOR b = fr TO tr 
430 IF z$(3,2) = “C” THEN LET 
v(2)=v(2) + 1: LET 
v(4) =v(4) 4- (v(3) < > 26) 

440 IF z$(3,2) = “R” THEN LET 
v(3) = v(3) + (v(3) < > 26): LET 
v(1)=v(1) + 1 

450 IF v(1) <25 AND v(2) <31 AND 
v(3) = 26 THEN GOTO 470 
460 IF v(1) >24 OR v(2) > 30 OR v(3) > 24 
OR v(4) >30 THEN GOTO 570 
470 IF v(1 ) < 1 OR v(2) <1 OR v(3) <1 OR 
v(4) < 1 THEN GOTO 570 
480 LET a$ = CHR$ (v(1 ) + 64) + STR$ 
v(2) + CHR$ (v(3) + 64) + STR$ v(4) 

+ o$ 

485 LET c = LEN a$: IFC > 8THEN RETURN 
490 RESTORE 1630: FOR q = 1 TO 11: LET 
f = 0: READ m$: F0Rw = 1 TO c 
500 IF m$(w) = “A” THEN GOSUB 1650: IF f 
THEN GOTO 560 

510 IF m$(w) = “N” THEN GOSUB 1670: IF f 
THEN GOTO 560 

520 IF m$(w) = “Z” THEN GOSUB 1710: IF f 
THEN GOTO 560 

530 IF m$(w) = “0” THEN GOSUB 1690: IF f 
THEN GOTO 560 

540 NEXT w: LET z = q: GOSUB 1140: 

IF NOT f THEN LET $$ = “□□□ 

□ □□□□”: FOR w = 1 TO c: LET 
s$(w) = a$(w): NEXT w: LET d$(b,a,9 TO 
16) = s$: LET d$(b,a,18) = CHR$ z: LET 
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d$(b,a,17) = “1”: NEXT b: NEXT a: 
RETURN 
550 GOTO 570 
560 NEXT q 
570 RETURN 

580 LET e = c: LET a$ = “DDD” 

590 PRINT # 1 ;AT 0,x; BRIGHT 1;“d” 

600 PAUSE 0: LET i = CODE INKEYS 
610 IF i> 88 THEN GOTO 600 
620 IF i = 13 THEN GOTO 650 
630 IF i = 12 THEN LET a$(4 — e) = 

LET e = e 4- 1 : LET x = x — 1 
640 LET a$(4 — e) = CHR$ i: PRINT #1;AT 
0,x;CHR$ i: LET x = x + 1 : LET e = e — 1 : 

IF e>0 THEN PAUSE 10: GOTO 590 
650 IF e>1 AND (d = 1 OR d = 4 OR d = 5) 
THEN GOTO 590 

655 IF e>0 AND (d = 2 OR D = 3) THEN 
GOTO 590 

660 PAUSE 10: PRINT #1;AT0,0; 

“□□□□□□□□□□□□□a 

□□□□□□□□□□□□□□□ 

□ □”:RETURN 
670 LET i$ = 

680 FOR z = 1 TO 3 
690 LET i$ = i$ + (a$(z) AND 
a$(z) <>“□”) 

700 NEXT z 

710 IF LEN i$ = 3 THEN IF i$(1)<“A” OR 
i$(1 ) > “X” OR i$(2)<“0” OR 
i$(2) >“9” OR i$(3) <“0” OR 
i$(3) >“9” THEN LET f = 1 : RETURN 
720 IF LEN i$ = 3 THEN IF VAL i$(2 TO 3) = 0 
ORVAL i$(2 TO 3) >30 THEN LET f=1: 
RETURN 

730 IF LEN i$ = 2 THEN IF i$(1 ) < “A” OR 
i$(1)>“X” OR i$(2) <“1” OR 
i$(2) >“9” THEN LET f=1: RETURN 
740 IF d = 2 THEN IF i$(1)< >“A” AND 
i$(1) < >“R” THEN LET f = 1 : RETURN 
750 IF d = 3 THEN IF i$(1)< >“C” AND 
i$(1)< >“R” THEN LET f = 1 : RETURN 
760 LET z$(d,2 TO ) = i$: LET z$(d,1) = CHR$ 
(LEN i$ + 48): LET f = 0: RETURN 
770 LET fc = (CODE z$(4,2)) — 64: LET 
tc = (CODE z$(5,2)) - 64: LET fr = VAL 
z$(4,3 TO (1 +VAL z$(4,1 ))): LET tr=VAL 
z$(5,3 TO (1+VAL z$(5,1 ))) 

780 IF z$(3,2) = “C” THEN IF fc< >tc OR 


fr > tr THEN LET f = 1 : RETURN 
790 IF z$(3,2) = “R” THEN IFfotcOR 
fr< >tr THEN LET f=1: RETURN 
800 LET f = 0: RETURN 
810 FOR y = 1 TO 30: FOR x = 1 TO 24: LET 
os = 0 

820 IF d$(y,x,17) = “1” THEN LET z = CODE 
d$(y,x,18): GOSUB 880: GOSUB 1010: 

LET st = 1 : LET a$ = STR$ t: LET os = LEN 
a$: IF t> 99999.99 THEN LET 
d$(y,x,17) = “5”: RETURN 
830 IF os>8 THEN LET st=os — 7 
840 IF os = 0 THEN GOTO 860 
850 LET s$ = “D □□□□□□□”: FOR 
u = st TO os: LET s$(u — st-1- 1 ) = a$(u): • 
NEXT u: GOSUB 1410: LET d$(y,x, TO 
8) = s$ 

860 LET i = IN 32766: IF i = 252 THEN PRINT 
#1;AT 0,0; PAPER 2; INK 7;“CALCULATE 
ABORTED”: RETURN 
870 NEXT x: NEXT y: RETURN 
880 LET s$ = d$(y,x,9 TO 16) 

890 IF z = 1 THEN LET v(1) = (CODE 
s$(1 )) — 64: LET v(2) =VAL s$(2): LET 
v(3) = (CODE s$(3)) — 64: LETv(4)=VAL 
s$(4): LET o$ = s$(5): RETURN 
900 IF z = 2 THEN LET v(1 ) = (CODE 
s$(1 )) - 64: LET v(2) = VAL s$(2 TO 3): 
LET v(3) = (CODE s$(4)) — 64: LET 
v(4) = VAL s$(5): LET o$ = s$(6): RETURN 

BE 

860 A$ = A$ + “D” 

870 AS = LEFT$(A$,LEN(A$) — 1 ) 

880 IF RIGHT$(A$,1) = “□” THEN 870 
890 PS = 1 :LN = LEN(A$):K = 0 
900 IF LN = 1 AND AS = “R” THEN FI = 3: 
RETURN 

910 IF LN = 1 THEN FI =0:RETURN 
920 IF LN = 0 AND RIGHT$(D$(R,C),8) = 
“□□□□□□□□” THEN A$ = 
CHR$(128):F1 =3:RETURN 
930 GOSUB790 
940 PS = PS + 1 + D1 
950 IF D1 =0 THEN FI =0:RETURN 
960 IF PS < LN — 1 THEN 930 
970 IF PS>LN THEN FI =3:RETURN 
980 Q1 $ = OP$:Q2$ = MID$(A$,PS,1 ): 
GOSUB 2670 


990 IF Q = 0 THEN FI =0:RETURN 
1000 D = VAL(RIGHT$(A$,1 )) 

1010 IF LN = PS THEN FI =3:RETURN 
1020 IF LN = PS + 1 AND D>1 AND D < 10 
THEN FI =3:RETURN 
1030 FI =0:RETURN 
1040 IF AS = “’’THEN RETURN 
1050 GOSUB 860 

1060 IF FI =0 THEN A$ = “”:INPUT 
“□ ENTER U ”;A$:G0T01 040 
1070 RETURN 
1080 PRINT“M”SPC(16) 

11 ,3 WORKING H” 

1090 FOR C = 1 TO CM 
1100 FOR R = 1 TO RM 
1110 D = 0:F = 0:OP = 8 
1120 F$ = LEFT$(D$(R,C),8) 

1130 Cl =ASC(MID$(F$,1,1)) — 64 
1140 IF Cl =64 OR Cl = -32 THEN 1230 
1150 IF Cl =18 THEN GOSUB 1320: 

GOTO 1230 

1160 R1 = VAL(MID$(F$,2,2)):IF R1 >9 THEN 
F = 1 

1170 C2 = ASC(MID$(F$,3 + F)) — 64: 

IF C2= — 32 THEN GOSUB 1360: 

GOTO 1230 

1180 IF C2= —3 THEN D = VAL(MID$ 

(F$,4 + F,2)):GOSUB 1360:GOTO 1230 
1190 R2 = VAL(MID$(F$,4 + F,2)): 

IF R2>9 THEN F = F + 1 
1200 Q1$ = OP$:Q2$ = MID$(F$,5 + F,1 ): 
GOSUB 2670 

1210 D = VAL(MID$(F$,6 + F,2)) 

1220 GOSUB 1360 
1230 NEXT R,C 
1240 RETURN 

1250 GET AS: IF A$ = “” THEN 1250 
1 260 IF AS < AA$ OR AS > BBS THEN 1 250 
1270 RETURN 

1280 GET A$:IF A$ = “” THEN 1280 
1290 IF A$ = AA$ THEN RETURN 
1300 IF AS = BBS THEN RETURN 
1310 GOTO 1280 
1320 AA$ = RIGHT$(D$(R, C),8) 

1330 IF RIGHT$(AA$,1) = “D” THEN 
AA$ = LEFT$(F$,7):G0T0 1330 
1 340 D$(R,C) = LEFT$(D$(R,C),8) + AA$ 
1350 RETURN 

1360 VI =VAL(RIGHT$(D$(R1,C1),8)):V2 = 

/ 







tF 
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VAL(RIGHT$(D$(R2,C2),8)) 

1370 ON Q GOSUB 1420,1430,1440,1450, 
1460,1470,1510,1550 
1380 AA$ = STR$(AN):IF LEN(AA$)>8 
THEN AA$ = “TOOL) BIG” 

1390 IF LEN(AA$) <8 THEN AA$ = “D” + 
AA$:GOTO 1390 
1400 D$(R,C) = F$ + AA$ 

1410 RETURN 

1420 AN = V1 +V2:RETURN 
1430 AN = V1 — V2:RETURN 
1440 AN = VI *V2: RETURN 
1450 AN = V1/V2:RETURN 
1460 AN = V2'V1/1 00:RETURN 
1470 AN = 0:FOR NN = R1 TO R2 
1480 AN = AN + VAL(MID$(D$ 
(NN,C1),9,8)) 

1490 NEXT NN 
1500 RETURN 

1510 AN = 0:FOR NN = C1 TO C2 
1520 AN = AN + VAL(MID$(D$(R1, 
NN),9,8)) 

1530 NEXT NN 
1540 RETURN 
1550 AN = V1:RETURN 
1560 PRINT“DO YOU WISH TO SAVE THIS 
DATA (Y/N)?” 

1570 AA$ = “Y”:BB$ = “N”: 

GOSUB1280 

1580 IF A$ = “N” THEN RETURN 
1590 INPUT"nNAME”;NM$: 
OPEN1,1,1,NM$ 

1600 PRINT# 1,RM:PRINT#1, CM: 

FOR R = 1 TO RM 
1610 FOR C = 1 TO CM 
1620 PRINT#1,CHR$(34) + D$(R,C) + 
CHR$(34) 

1630 NEXT C,R:CLOSE1:RETURN 
1640 PRINT u QH|JDO YOU WISH TO 
LOAD” 

1650 PRINT“||AN EXISTING FILE (Y/N)?” 
1660 AA$ = “Y”:BB$ = “N”:GOSUB1280 


a 

870 DEF FNformat 

880 IF LENA$>8 A$ = LEFT$(A$,8) 

890 IF Type = 8 AND LENA$ < 8 REPEAT 
A$ = A$ + “ □ ”:UNTIL LENA$ = 8 
900 IF Type = 0 AND LENA$<8 REPEAT 
A$ = “D” + A$:UNTIL LENA$ = 8 
910 =A$ 

920 DEF FNboxcheck(pos) 

930 LOCAL a$,b$,p,d 
940 a$ = MID$(A$,pos,3) 

950 b$ = LEFT$(a$,1 ) 

960 IF b$<“A” OR b$>“X” =0 
970 p = VAL(RIGHT$(a$,2)) 

980 d = 2:IF p < 1 0 d = 1 

990 IF p > Cols OR p<1 =0 

N.B. In Line 990, Electron users should 

change Cols to ROWS 

1000 =d 

1010 DEF FNCheck 

1020 LOCALpos,ln,f,d,k 

1030 A$ = A$ + “D” 

1040 REPEAT A$ = LEFT$(A$,LENA$ — 1) 
1050 UNTIL RIGHT$(A$,1) < >“□” 

1060 pos = 1:ln = LENA$ 

1070 IF In = 1 AND A$ = “R” =3 
1080 IF In = 1 =0 

1090 IF ln = 0 AND RIGHT$(D$(Row,Col), 

8) = STRING$(8,‘‘D”) A$ = CHR$(128): 

= 3 

1100 REPEAT 
1110 k = k + 1 
1120 f=FNboxcheck(pos) 

1 1 30 pos = pos + 1 + f 

1140 UNTIL f = 0 OR pos> = In OR k = 2 

1150 IF f=0 =0 

1160 IF pos > In =3 

1170 IF INSTR(Op$,MID$(A$,pos,1 )) = 0 =0 
1180 d = VAL(RIGHT$(A$,1 )) 

1190 IF In = pos =3 

1200 IF ln = pos + 1 AND d> — 1 AND 


1260 f = FNCheck 

1270 IF f=0 PRINTTAB(0,vpos)CHR$129; 
‘‘INCORRECT:PLEASE RE-DOD”;:INPUT 
“”A$ 

1280 UNTIL f>0 

1290 ENDPROC 

1300 DEF PROCcalculate 

131 0 PRINTCHRS1 29;CHR$1 36;“WORKING” 

1320 FOR c% = 1 TO Cols 

1330 FOR r% = 1 TO Rows 

1340 d% = 0: f% = 0:op% = 8 

1 350 F$ = LEFT$(D$(r%,c%),8) 

1360 c1% = FNc(1) 

1370 IF cl % = 64 OR c1%= -32 GOTO 1470 
1400 r1% = FNr(2):IF r1%>9 f% = 1 
1410 c2% = FNc(3 + f%):IF c2%= -32 
PR0Carith:G0T01 470 
1420 IF c2%= —3 d% = FNr(4 + f%): 

PROCarithrGOTOI 470 
1430 r2%= FNr(4 + f%):IF r2%>9 
f% = f% + 1 

1440 op%= INSTR(Op$,MID$(F$,5 + f%,1 )) 
1450 d%= FNr(6 + f%) 

1460 IF op% < >0 THEN PROCarith 
1470 NEXT 
1480 NEXT 
1490 ENDPROC 

1500 DEF FNr(p%) =VAL(MID$ 

(F$,p%,2)) 

1510 DEF FNc(p%) = ASC(MID$ 

(F$,p%,1))— 64 

1520 DEF FNvl =VAL(RIGHT$(D$(r1%, 
c1%),8)) 

1530 DEF FNv2 = VAL(RIGHT$(D$(r2%, 
c2%),8)) 

1540 DEF P R 0 Cgetbet ( a$ , b$ ) 

1550 REPEAT A$ = GETS 
1 560 UNTIL A$ > a$ AND A$ < b$ 

1570 ENDPROC 

1580 DEF PROCgetabs(a$,b$) 



1670 IFA$ = ‘‘N” THEN RETURN 
1680 INPUT“NAME”;NM$:OPEN 1,1, 
0,NM$ 

1690 INPUT#1,RM,CM:FOR R = 1 TO 
RM:FOR C = 1 TO CM 
1700 INPUT# 1,D$(R,C):NEXTC,R 


d < 1 0 =3 
1210 =0 

1220 DEF PROCformulacheck 
1230 LOCALvpos,f 

1 240 vpos = VPOS + 1 : 1 F AS = TH EN 

ENDPROC 
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1590 REPEAT A$ = GET$:UNTIL A$ = a$ OR 
A$ = b$ 

1600 ENDPROC 
1670 DEF PROCarith 
1680 LOCALans,a$,d 

1690 ON op% GOSUB 1780,1790,1800,1810, 
1820,1830,1870,1910 
1700 IF d% = 0 THEN ans = INT(ans + .5): 
GOTO 1720 

1 71 0 @% = &1 020008 + d%*256 
1720 a$ = STR$(ans):IF LENa$ > 8 
a$ = CHR$135 + “TOO BIG” 

1730 IF LENa$<8 a$ = “(H” + a$: 
GOTO1730 

1740 D$(r%,c%) = F$ + a$ 

1750 @% = &90A 
1760 ENDPROC 

ma 

660 A$ = D$(I,J):B$ = MID$(A$,2) 

670 AT=ASC(A$) 

680 IF AT = 1 28 THEN PRINTSTRING$(7,32);: 
GOTO 720 

690 IF AT = 129 OR AT=130 THEN PRINT 
USING “%□ □ □ □ □%”;B$;: 

GOTO 720 

700 FOR U = 1 TO LEN(B$):IF MID$ 

(B$,U,1) < >CHR$(32) THEN PRINT 
MID$(B$,U,1); 

710 NEXTU 
720 RETURN 

730 Cl = ASC(Z$) — 64:C2 = VAL(MID$ 
(Z$,2)):V = D(C1 ,C2):RETURN 
740 PRINT@448, “WORKING” 

750 FOR J = 1 TO RX 
760 FOR 1 = 1 TO CX 

770 D(I,J) = 0:IF ASC(D$(I,J)) = 129 THEN 
D(I,J) = VAL(MID$(D$(I,J),2» 

780 NEXT l,J 
790 FOR J = 1 TO RX 
800 FOR 1 = 1 TO CX 

810 PRINT@448, “WORKING ON CELLD”; 

CHR$(I + 64);MID$(STR$(J),2) 

820 IF ASC(D$(I,J)) < > 131 THEN 1130 
830 A$ = MID$(D$(I,J),2) 

840 0$ = MID$(A$,7,1) 

850 IF 0$ = “&” THEN 1050 
860 IF 0$ = “$” THEN 1090 
870 Z$ = LEFT$(A$,3) 


880 GOSUB 730 

890 VI =V:Z$ = MID$(A$,4,3):G0SUB 730: 
V2 = V 

900 DP = VAL(RIGHT$(A$,1 )) 

910 ON INSTR(1,0P$,0$) GOSUB 
1000,1010,1020,1030,1040 
920 OV = 0:IF DP = 0 THEN PU$ = 
“#######”: MP = 7: 

GOTO 950 

930 PU$ = STRING$(7 — (DP + 1 ),“#”) + 
“.” + STRING$(DP,“#”):MP = 7 — 

(DP + 1) 

940 IF LEN(PU$) >7 THEN RV$ = 

“□ □ <0V>”:0V = 1 
950 D(I,J) = RV 

960 IF RV < 0 THEN MP = MP-1 
970 ML = LEN{MID$(STR$(INT(RV + .5)), 2)) 
980 IF ML>MP THEN RV$ = “DD 
<0V>”:0V = 1 
990 GOTO 1160 
1000 RV = VI +V2:RETURN 
1010 RV = VI — V2:RETURN 
1020 RV = VI *V2:RETURN 
1030 IF V2 = 0 THEN RV = 0:RETURN ELSE 
RV = V1/V2:RETURN 
1040 RV = VI 'V2/1 00:RETURN 
1050 PI = ASC(A$) — 64: P2 = ASC(MID$ 
(A$,4,1 )) - 64:C2 = VAL(MID$(A$,2,2)): 
RV = 0 

1060 FOR Cl = PI TO P2 
1070 RV = RV + D(C1,C2):NEXT 
1080 DP = VAL(RIGHT$(A$,1 )): 

GOTO 920 

1090 PI = VAL(MI D$(A$,2,2)):P2 = VAL 
(MID$(A$,5,2)):C1 =ASC(A$) 

— 64:RV = 0 

1100 FOR C2 = P1 TO P2 
1110 RV = RV + D(C1,C2):NEXT 
1120 DP = VAL(RIGHT$(A$,1 )): 

GOTO 920 

1130 IF ASC(D$(I,J)) < >128 
THEN 1150 

1140 RV$ = STRING$(7,32):GOTO 1160 
1150 RV$ = MID$(D$(I,J),2) 

1160 IF I > =CS AND l< = CS + 3 AND 
J> = RS AND J<=RS + 11 THEN 
PRINT@(J — RS)'32 
+ 35 + (I — CS)'7,“”;:PF = 1 
ELSE PF = 0 


1170 IF (ASC(D$(I,J)) >= 128 AND 
ASC(D$(I,J)) < =130) OR 0V=1 THEN 
1200 

1180 IF PF = 1 THEN PRINT USING 
PU$;RV; 

1190 GOTO 1210 

1200 IF PF = 1 THEN PRINTUSING 
“%□ □ □ □ □%”;RV$; 

1210 NEXT I, J 
1220 RETURN 

1230 CLS:INPUT “DO YOU WANT TO SAVE 
THIS SHEETD □(Y/N)”;A$ 

1240 IF A$o“Y” THEN 1340 
1250 LINE INPUT “FILENAME:”;F$ 

1260 OPEN “0”,# — 1,F$ 

J 270 FOR J = 1 TO RX 
1280 FOR 1 = 1 TO CX 
1290 IF ASC(D$(I,J)) = 128 THEN 1320 
1300 Z$ = D$(I,J):MID$(Z$,1,1 ) = CHR$ 
(ASC(MID$(Z$,1 ,1 )) — 95) 

1310 PRINT# -1,STR$(I),STR$(J): 

PRINT# -1,Z$ 

1320 NEXTI,J 
1330 CLOSE# -1 

1340 CLS:M0 = 1:G0SUB 70:RETURN 
1350 CLS:INPUT “DO YOU WANT TO LOAD A 
SHEET FROMTAPE (NOTE SHEET IN 
MEMORY WILL BE MERGED WITH THAT 
ON TAPE) □ □ □ □ (Y/N)”;A$ 

1360 IF A$< >“Y” THEN 1480 
1370 PRINT“PRESS enter TO LOAD NEXT FILE 
ONTAPE- OR INPUT FILENAME 
NOW”:PRINT 

1380 LINE INPUT “FILENAME:”;F$ 

1390 OPEN “I”, #— 1,F$ 

1400 IF EOF(-I) THEN 1470 
1410 INPUT# — 1,A$,B$:LINE 
INPUT# -1,C$ 

1420 MID$(C$,1,1) = CHR$(ASC(MID$ 
(C$,1,1)) + 95) 

1430 Cl = VAL(A$):C2 = VAL(B$):D$(C1 , 

C2) = C$ 

1440 IF Cl >CX THEN CX = C1 
1450 IF C2>RX THEN RX = C2 
1460 GOTO 1400 
1470 CLOSE# -1 

1480 CLS:CC = 1 :CR = 1 :CS = 1 :RS = 1 : 

M0 = 1:G0SUB 70:RETURN 
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CUMULATIVE INDEX 


An interim index will be published each week. There will be a complete index in the last issue of INPUT. 


A 

Animation 

of UDGs in cliffhanger 992-997 

using colour fill techniques 
Acorn 955-959 

using GCOL 3 

Acorn 999-1000 

using paged graphics 

1022-1027, 1132-1137 

Applications 

calendar and diary program 

1010-1016, 1017-1021, 1064-1067 
hobbies file, extra options 947-952 

magnification program 1081-1087 

spreadsheet program 

1118-1126, 1172-1176 
text-editor program 

852-856, 878-883, 914-920 


B 

BASIC 

adding instructions to 

Acorn, Dragon, Spectrum 844-851 
Basic programming 
analyzing and storing sounds 1091-1095 
animation with paged graphics 


1022-1027, 1132-1137 
colour commands, Acorn 953-959 

Computer Aided Design 998-1004 
designing a new typeface 838-843 

drawing conic sections 859-863, 889-895 
how programs are stored 1 106-1 1 12 
mathematics of growth 1049-1056 

mechanics, principles of 933-939 

multi-key control 974-979 

musical chords and harmonies 985-991 
patterns from nature 1 164-1 171 

prediction by computer 1 1 58-1 1 63 
programming function keys 825-829 
secret codes 960-965,1044-1048 

sound envelopes 

Acorn, Commodore 64 1138-1144 

speeding up BASIC programs 921-927 


c 

Calendar and diary program 

1010-1016, 1017-1021, 1064-1067 

Chords, musical 

definition 985-986 

programs to play 

Acorn, Commodore 64 986-991 

Cliffhanger game 

part 1— title page 904-913 

part 2— adding instructions 928-932 

part 3 — adding a tune 966-973 

part 4 — graphics and merging 992-997 

part 5 — setting the scene 1034-1043 

part 6 — perils and rewards 1057-1063 

part 7— initializing routine 1 101-1 105 

part 8 — synchronizing routine 

1127-1131 

part 9 — scoring routine 1 145-1151 

Codes, secret 960-965, 1044-1048 

Colour 

defining in machine code 1034-1043 


filling in with 

Acorn 953-959 

in Teletext mode 

BBC 1068-1073 

routines for changing 

Commodore 64 872-877 

Computer Aided Design 
rubber-banding and picking 
and dragging 998-1004 

Conic sections 857-863, 889-895 


D 

Digital clock routine 896-898 

Dynamics, programs to 
illustrate 1164-1171 


E 

Envelope, 

of orbits 1164-1167 

sound 

Acorn, Commodore 64 

968-971, 1138-1144 
in musical harmony programs 986-991 


F 

Filling in with colour 

Acorn 953-959 

Fish population program 1 1 70-1 1 7 1 

Fox and geese game 
part 1 — principles and graphics 

1096-1100 

part 2 — initializing and 

mapping the moves 1113-1117 

part 3— higher levels 1 152-1157 

Fruit machine game 

1028-1033, 1074-1080 
Function keys, programming 
Acorn, Commodore 64, Vic 20 826-829 


G 

Games 

cliffhanger 904-913, 928-932, 966-973, 
992-997, 1034-1043, 1057-1063, 1101-1105, 
1127-1131, 1145-1151 

fox and geese 

1096-1100, 1113-1117, 1152-1157 
fruit machine 1028-1033, 1074-1080 
goldmine 830-837, 864-871 

lunar touchdown 1088-1090 

magnification 1081-1087 

multi-key control for 974-979 

othello 980-984, 1005-1009 

wordgame 899-903, 940-945 

Goldmine game 830-837, 864-871 

Graphics 

colour commands, Acorn 953-959 

effects using curves 857-863, 889-895 
hi-res 

for custom typeface 838-843 

setting up new commands 

Commodore 64 872-877 

magnification program for 1081-1087 


paged, for animation 

1022-1027, 1132-1137 


patterns from nature 1 164-1 171 

picking and dragging 1000-1004 

rubber-banding 998-1000 

trace of sound 1092-1095 

using Teletext mode, BBC 1068-1073 
Growth, measuring 1049-1056 


H 

Hobbies file, extra options for 947-952 

I 

Instructions, adding to BASIC 

Acorn, Dragon, Spectrum 844-851 


K 

Keypresses, multiple, programming for 

974-979 


L 

Letter-generator program 838-843 
Lunar touchdown game 1088-1090 

M 

Machine code 

games programming 
see cliffhanger 

merging routines 992-997 

routines for hi-res graphics 

Commodore 64 872-877 

routine to alter BASIC 844-849 

timer routine 896-898 

tune routine 966-973 

Magnification program 1081-1087 

Mathematical functions 

in mechanics 935 

in spreadsheet program 1120 

speedy use of 923-924 

to assess population tendencies 

1170-1171 

to draw curves 857-863, 889-895 

to draw patterns from orbits 

1164-1170 

to measure growth 1049-1056 

Mechanics 

programs to show principles of 933-939 

Memory 

how BASIC programs are stored in 

1106-1112 

paged graphics in 1023-1027, 1132-1137 
Multi-key control, programming for 

974-979 

Music 

analyzing and storing 1091-1095 

chords and harmonies 985-991 

machine code routine for 966-973 

N 

Numbers 

Fibonacci 1056 


generation program 

1054-1056 

0 


Orbits, patterns from 

1164-1171 

Othello board game 


980-984, 

1005-1009 

Q_ 


Paged graphics 1023-1027, 

1132-1137 

Patterns from nature 

1164-1171 

PLOT 


new commands, Acorn 

953-959 

Pools simulation program 

1158-1160 

Prediction by computer 

1158-1163 

R 


RND function 


in computing probability 

1160-1163 

Robotics 

884-888 

s 


Search routines 


binary and serial 

924-927 

in text-editor program 

914-920 

single pass 

1162-1163 

Sounds 


analyzing and storing 

1091-1095 

envelopes for modifying 


Acorn, Commodore 64 

1138-1144 

Sort routines 


in hobbies file program 

947-953 

in text-editor program 

914-920 

Speeding up BASIC programs 

i 921-927 

Spreadsheet program 


part 1 

1118-1126 

part 2 

1172-1176 

T 


Teletext mode, BBC 

1068-1073 

Text-editor program 


part 1 — basic routines 

852-856 

part 2 — editing facilities 

878-883 

part 3 — sorting, searching, 


formatting and printout 

914-920 

Timer routine 


for BASIC lines 

92l2 

machine code 

896-898 

Typeface, setting up new 

838-843 


V 

Variables 

managing for program speed 923-925 
setting in machine code game 1 1 27-1 1 3 1 
storing in memory 1 1 06-1 1 1 2 

w 

Waveforms 

displaying and storing 1092-1095 

modulation of 1 1 38-1 1 39, 1 142 

Wordgame 899-903, 940-945 


The publishers accept no responsibility for unsolicited material sent for publication in INPUT. All tapes and 
written material should be accompanied by a stamped, self-addressed envelope. 




COMING IN ISSUE 38. . . 


LJ A fearsome Martian spider, a neurotic 
window cleaner, a bevy of balloons and 
an accumulation of arrows are all part of 
FREDDY AND THE SPIDERS FROM 
MARS, our BASIC arcade game 

U Complete the SPREADSHEET 
program and plan your finances 

Li Got ten thumbs? Can’t draw a thing? 
GRAPHICS OF ROTATION allows you 
to create realistic three-dimensional 
views of any shape you choose 

O Use statistical techniques to MODEL 
REALITY. Find out how to apply 
statistical models 

— /Sight some circling seagulls in 
CLIFFHANGER on the SPECTRUM, 
COMMODORE and ACORN machines 
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